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ABSTRACT 

A  two-dimensional,  incompressible  model  is  presented  for  the  response 
of  laminar  and  turbulent  attached  jets  to  arbitrary  control  pressures 
applied  to  either  or  both  control  ports  of  a  flueric  bistable  amplifier. 
Effects  of  the  nonlinear,  lumped-parameter,  inductive  and  resistive  sup¬ 
ply,  control,  vent,  and  output  lines,  and  effects  of  the  wall  opposite 
the  attachment  wall  are  considered.  All  geometric  parameters  such  as 
wall  angle  and  offset,  vent  location,  line  size,  and  nozzle  width  may 
be  varied. 

Switching  time  data  for  seven  different  geometries  and  steady-state 
attachment-point  data  for  both  closed  and  open  control  ports,  as  well 
as  the  response  of  an  actual  amplifier  compare  favorably  with  analytical 
predictions.  The  response  of  a  laminar  jet  to  a  ramp  input  pressure 
change  is  also  given.  The  analytical  results  are  in  good  agreement  with 
experimental  data.  These  results  indicate  that  the  analytical  model  is 
valid  and  that  it  can  be  used  as  a  basis  in  designing  future  flueric 
amplifiers , 
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1.  INTRODUCTION 


A  fundamental  component  of  many  fluidic  systems  is  the  wall-attach¬ 
ment  amplifier.  In  spite  nf  its  importance  and  the  great  deal  of  work 
expended  on  it,  the  design  of  this  amplifier  has  been  based  primarily 
on  cut-and-try  techniques  because  a  sufficiently  accurate  analytical 
model  has  not  been  available.  The  advance  toward  a  design  theory  of  a 
wall-attachment  amplifier  began  with  Bourq\ie  and  Newman^  who  used  the 
attachment  bubble  control  volume  concept  to  determine  the  point  of  re¬ 
attachment.  T.iis  concept,  with  many  variations,  has  proved  fruitful, 
and  by  its  use,  the  theoretical  design  has  been  extended  by  Wada  and 
Shimizu^  and  Kitiiura  cind  Mitsuoka.’  Time-dependent  analyses  have  been 
presented  recently  by  Epstein,"*  Wilson, '  Lush,®“®  and  Goto  and 
Drzewiecki.’  The  re-'^erences  cited  by  Epstein  adequately  summarize  the 
previous  work . 

In  this  paper,  we  present  a  practical  analytical  model  for  the 
internal  dynamics  and  switching  characteristics  of  a  flueric  wall- 
attachment  amplifier  with  a  sharp  splitter.  This  model  improves  on 
our  previous  one,®  because  the  affects  of  the  splitter  and  outputs  are 
considered  as  well  as  the  fact  that  the  jet  is  allowed  to  attach  to  the 
opposite  wall.  In  addition  to  discussing  the  improved  model,  we  also 
give  an  account  of  our  previous  models.  We  believe  that  thi.s  is  worth¬ 
while  for  two  reasons:  (1)  the  steps  through  which  the  present  model 
has  evolved  are  clarified,  and  (2)  the  theory  and  data  on  the  simple 
model  should  prove  useful  in  formulating  theories  for  other  fluidic 
designs . 

The  analysis  for  the  inner  portion  of  the  amplifier  is  based  on  a 
fully  developed  two-dimensronal ,  incompressible  laminar  or  turbulent 
jet.  It  is  reasonable  to  postulate  fully  developed  jet  flow  in  the 


^  Bourque,  C.  and  Newman,  B.G.,  "Reatcaohment  of  a  Two-Dimensional,  Incompressible 
Jet  to  an  Adjacent  Flat  Plate,"  The  Aeronautical  Quarterly,  Vol.  XI,  August  1960, 

^  Wada,  T.  and  Shimizu,  A.,  "Experimental  Study  of  Attaching  Jet  Flow  on  an  Inclined 
Flat  Plate  with  Small  Offset,"  Proc.  2nd  IFAC  Symposium  on  Fluidics,  Prague,  1971. 

®  Klmura,  M.  and  Mitsuoka,  T.,  ".Analysis  and  Design  of  Wall  Attachment  Devices  by  a 
Jet  Model  of  Unsymmetrical  Velocity  Profile,"  Proc.  IFAC  Symposium  on  Fluidics, 
London,  1968. 

^  Eostein,  M. ,  "Theoretical  Investigation  of  the  Switching  Mechanism  in  a  Bistable 
Wall  Attachment  Fluid  Amplifier,"  ASME  Publication  70-F1CS-3,  June  1970,  also  J. 
Basic  Eng.,  pp .  55-62,  March  1971. 

*  Wilson,  M.P.,  Jr.,  "The  Switching  Process  In  Bistable  Fluid  Amplifiers,"  ASME 
Publica  ion  69-Flcs-28,  June  1969. 

*  Lush,  P.A.,  "Investigation  of  the  Switching  Mechanism  in  a  Large  Scale  Model  of 
a  Turbulent  Reattachment  Amplifier, "  Second  Cranfleld  Fluidics  Conference,  Cam- 
Bridge,  England,  1967. 

’’  Lush,  v.A.  ,  "A  Theoretical  and  Experimental  Investigation  of  the  Switching  Mechan¬ 
ism  in  a  Wall  Attachment  Fluid  Amplifier,"  Proc.  IFAC  Symposium  on  Fluidics, 
London,  1968. 

®  Lush,  P.A.,  "The  Development  of  a  Theoretical  Model  for  the  Switching  Mechanism  of 
a  Wall  Attachment  Fluid  Amplifier,"  Ph.D.  Thesis,  University  of  Bristol,  Dept,  of 
Aero.  Engr.,  September  1963. 

®  Goto,  J.M.  and  Drzewleckl,  T.M.,  "Reattached  Jet  Response  to  Input  Pressure  in  a 
Non-Loaded  Fluidic  Bistable  Configuration,"  Fifth  Cranfleld  Fluidics  Cotiference, 
Uppsala,  Sweden,  June  1972, 
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interaction,  chamber  as  long  as  control  port  gaps  are  present.  Note, 
however,  that  prior  analyses  assuming  a  fully  developed  laminar  jet 
profile  with  no  control  port  gap  nave  grossly  overestimated  the  attacli- 
ment  point  distance  for  geometries  having  small  wall  angles  and  zero 
offset  (Comparin,  et  al.'^).  When  there  is  no  control  port  gap,  one 
might  think  that  the  flow  never  fully  separates  and  therefore  attaches 
very  quickly.  Qualitative  water  table  observations  show  that  when  a 
control  gap  is  provided,  the  attachment  point  moves  far  downstream, 
more  in  agreement  with  the  fully  developed  jet  predictions.  Thus, 
the  use  of  a  fully  developed  jet  analysis  is  justified  only  when  con¬ 
trol  port  gaps  exist.  Similar  arguments  apply  to  the  turbulent  jet. 

It  is  interesting  to  note  that  the  geometric  scaling  laws  for  turbulent 
confined  flows  are  unknown,  so  that  for  turbulent  jets,  the  analysis 
should  hold  only  when  fully  developed  turbulence  is  established.  As 
long  as  that  criterion  is  met,  or  when  the  jet  is  laminar,  there  should 
not  be  any  physical-size  limitations  on  the  devices  for  which  the  theory 
holds . 

The  outer  portion  of  the  amplifier  model  consists  of  channels  for 
the  supply,  control,  vent,  and  output  linos.  The  lines  are  character¬ 
ized  as  lumped-parameter  inductive  and  nonlinear  resistive  components. 
Tile  lines  are  coupled  to  tlie  inner  region  pressures;  therefore,  flows 
through  all  lines  vary  with  changes  of  jet  position.  Most  of  the  basic 
assumptions  have  been  used  extensively  by  previous  authors. Addi¬ 
tional  assumptions  about  the  nature  of  the  flow  were  necessary  for  this 
model  because  variable  line  flows  and  the  vortex  in  the  attachment  bub¬ 
ble  were  considered.  Furthermore,  low-offset  cases  may  be  studied  by 
allowing  ti.e  jet  to  flow  along  the  opposite  wall  as  opposeid  to  terminat¬ 
ing  the  proce''^  as  soon  as  the  jet  touches  the  opposite  wall."'®  A 
sharp  splitter  allows  the  assumption  of  a  lossless  peeling  off  of  raomeri- 
tum  flux  from  the  jit,  and  the  motion  of  the  jet  past  the  controls  per¬ 
mits  the  calculation  of  control  resistance  based  on  the  jet-wall  spacing, 
if  necessary. 

2 .  FORMULATION 

The  analysis  is  formulated  for  the  model  shown  rn  figures  1  and  2. 
Distances  x  and  y  are  referred  to  the  geometric  axis  of  syrrm;etry  and  the 
supply  nozzle  exit  plane,  whereas  n  and  s  are  referred  to  the  curving 
jet  centerline.  A  distance  n  is  measured  normal  to  the  jet  centerline 
at  some  jet  centerline  arc  length  s.  The  Goertler  velocity  profile 
(Schlichting '  p.  606)  is  assumed  for  the  turbulent  supply  jet  and  the 
Schlichting  profile  (Schlichting , ' *  p.  167)  for  the  laminar  jet.  The 
turbulent  profile  is  independent  of  Reynolds  number  (Re)  but  contains  a 
jet-spread  parameter,  \/hereas  the  laminar  profile  is  Reynolds-number-de- 
pendent,  containing  no  undetermined  parameters. 

The  centerline  of  the  jet  is  chosen  as  the  streamline  mo.st  signifi¬ 
cant  to  the  downstream  temporal  position  of  the  jet.  In  general,  the 
subscripts  a  and  i  denote  the  active  (attached)  and  inactive  (opposite) 
side  of  the  flueric  amplifier,  respectively.  The  extended  ]Gt  center¬ 
line  intersects  the  attachment  wall  at  the  point  denoted  by  the  sub¬ 
script  w.  All  unprimed  quantities  are  dimensionless.  Quantities  with 
dimensions  of  length  are  normalized  with  respect  to  the  supply  nozzle 
width  b+*  (The  subscript  +  refers  to  an  invariant  supply,  source,  or 


*  °C:ompai-in,  R.A.  ,  Jenkins,  W.C.  and  Moore,  K.B.,  ‘Met  Rent  tachraent  at  Low  Reyn  ’  's 
Number  and  Moderate  Aspect  Ratios,"  ASMK  Publication  67-Fi;-25,  May  1967. 

’'schlichting,  H.  ,  Boundary  Laver  Theory.  MoLraw-Hi i 1 ,  Now  York,  1960. 
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b  GENERAL  FLOW  MODEL 


Figure  I.  Attachment  model 


reference  condil.ion.)  Areas  and  volumes  are  per  unit  depth  and  are 
normalized  by  b|  and  b|^  respectively.  All  pressures  are  gage.  Static 
pressures  are  normalized  with  respect  to  the  supply  total  pressure  p+. 
Since  the  supply  nozzle  exit  pressure  pg  is  not  a  constant  in  this 
model/  the  supply  flow  per  unit  depth  Qg  is  not  constant,  and  so  for 
volumetric  flow  normalization  purposes,  an  arbitrary  constant  flow 
is  defined  ny  the  following  relationship  to  the  actual  supply  total 
pressure  pj.: 

q;  =  b;/2p;/p-  (i) 

where  p'  is  the  fluid  density.  Velocities  u'  are  normalized  with  re¬ 
spect  to  =  Q|/b^.  Note  that  flow  is  per  unit  depth.  The  Reynolds 
number  Re^  is  defined  as 


where  v  is  the  kinematic  viscosity.  Subsequent  equations  involving  I 
flow  and  velocity  will  differ  from  those  of  previous  analyses  by  some  ■ 
factor  of  the  supply  flow  Q|  .  Time  t'  is  normalized  with  =  t4-  ■ 
Ordinarily,  t^  is  the  transport  time  based  on  the  supply  nozzle  1 
width  (i.e.,  the  time  required  for  a  fluid  particle  moving  at  n 
the  supply  jet  exit  velocity  Ug  to  travel  a  distance  of  one  nozzle  || 
width).  However,  b^-ZUg  =  D^^/Qg  is  the  actual  transport  time,  but  it  I 
has  lost  its  significance  because  it  varies  in  a  time-dependent  I 
process.  a 


Figure  lb  shows  all  the  flows  and  their  respective  pressure  differ¬ 
ences,  The  bvibble  and  opposite-side  pressures  are  represented  by  mean 
pressures  p^,  and  p^,  and  compensation  for  the  bubble  pressure  distribu¬ 
tion  is  made  by  assuming  that  the  active  control  exit  pressure  is  equal 
to  Pav  =  (Pb  Pi)/2,  The  supply  jet  exits  to  an  interaction  region 
that  has  a  pressure  lower  than  ambient,  and  it  is  assumed  that  the  exit 
pressure  of  the  supply  can  be  represented  by  an  average  pressure 
Ps  “  (pav  +  Pi)/2  (see  fig.  lb).  Thus,  the  supply  and  control  flows 
are  coupled  to  the  interaction  region  by  these  pressures. 

2.1  Supply  Jet  Velocity  Profile  Equations, u 

In  terms  of  the  coordinates  s  and  n,  the  normalized  form  for  the 
fully  developed  two-dimensional  incompressible  jet  velocity  profile  is 

u  =  ki  sech^  k2n  (2) 


(3a) 


(3b) 


Cl  =  0,4543  Cz  =  0.2752 


wher^ 


hi 

kz 

ki 

k2 


=  Q 


3a 


^14 (s  +  s^q) 


Va 


s  +  s. 


VO 

'1/3 


+  ^vo) 


TURBULENT 


>  LAMINAR 


and  s  is  the  distance  measured  along  the  curving  jet  centerline.  The 
unknowns  in  eq  (2)  and  (3)  are  the  supply  flow  Qg,  the  point  source 
distance  3^^^,  and  the  fluid  velocity  u. 

The  turbulent  jet  spread  parameter  a  must  be  specified.  When  ap¬ 
plied  to  attached  jets,  it  is  usually  assigned  a  value  most  consistent 
with  experimental  data.  The  work  of  I.ush®“®  and  Kimura  and  Hitsuoka 
indicates  that  7.67<a<14.  A  constant  value  of  10  has  been  chosen 
for  use  in  this  model, 

2 . 2  Attachment  Bubble  VolU'.ie,  Vb 

The  jet  centerline  arc,  the  offset,  and  the  attacl.raent  walls  define 
the  contro  i^olume  of  the  attachment  bubble  (fig.  lb)  .  The  flows 
crossing  the  control  volume  boundaries  are  one-half  of  the  total  supply 
flow  Qg/2 ,  the  control  flow  Q^a'  khe  vent  flow  khe  downstream 

flow  Q downs tr earn*  There  is  no  vent  flow  Q^a  into  the  bubble  if  the  jet 
is  attached  upstream  of  the  vent.  The  downstream  flow  is  determined  by 
going  inside  the  control  voluiae  and  calculating  the  difference  between 
returned  flow  Qr  and  the  sum  of  the  half-supply  flow  Qs/2  and  entrained 
flow  Qe-  The  differential  equation  governing  the  time  rate  of  change  of 
the  bubble  volume  Vb  is  the  continuity  equation 


dt  2  ^  ^ca  '''  "^va  ^downstream 


1 1 


but 


SO  thct 


^downstream  2  ^ 

dVjj 

dt"  ~  ^ca  “  Qe  Qr  ■*■  ^va 


(4) 


This  is  then  a  capacitive  effect  caused  by  bubble  'olume  change. 


Tne  difference  between  the  returned  and  entrained  flows  is  obtained 
by  integrating  the  assumed  jet  velocity  profile  from  the  entrainment 
streamline  to  the  attachment  streamline  at  the  attachment  u-iint.  This 
results  in 


Qe 

^e 


“  ^r  2  ®vo)  1 

“  Qr  =  ^  ^ 
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(5a) 

(5b) 


where  s„  is  the  jet  centerline  arc  length  from  the  nozzle  exit  to  the 
wall  and  the  attachment  parameter  T„  is  defined  by  aq  (6). 
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(6a) 

(6b) 


where  ny  is  the  distance  from  the  jet  centerline  to  the  attachment 
streamline  evaluated  at  the  arc  distance  s^  and  is  a  function  of  the 
attachment  angle/  9  .  The  unknowns  in  eg  (5)  and  (6)  are  the  distances 
Syo#  and  s^. 


2 . 3  Jet  Deflection  Equation 


The  equation  for  the  jet  deflection  angle  6  is  obtained  by  combining 
the  continuity,  Bernoulli,  and  momentum  equations  in  the  interaction  re¬ 
gion  and  assuming  t'  at  the  supply  and  control  flows  interact  without 
losses.  Assuming  that  the  width  of  the  supply  jet  b4.  is  constant  in  the 
interaction  region,  then  from  continuity,  the  jet  vel.ocity  ug  is  con¬ 
stant,  and  by  Bernoulli's  equation  the  supply  jet  pressure  p^  throughout 
the  interaction  is  also  constant.  The  control  pressure  p^y  and  pi,  pro¬ 
ducing  forces  acting  normal  to  the  supply  jet,  and  the  control  momenta, 
p'u,lbc,  acting  in  the  y  direction,  cause  the  jet  to  deflect  through  an 
angle  B ,  as  shown  in  figure  lb.  The  force  balance  equation  in  the  lon¬ 
gitudinal  (x)  direction  is 


av 


pj)bjL  tan 


ps^;  - 


B  =  p'u^-b_j.  cos  6  -  p'u^’n]. 


Normalizing  with  p4.t>+  (-  P'Q4.^/2b|)  and  rearranging  gives 

<Pav  -  Pi)bctan  6  +  (pg  +  2Qg")  =  (pg  +  2QgMcos  6  {7a) 


By  introducing  the  transverse  forces  Fg  and  Ft ,  made  up  by  summing  the 
control  momenta  2QcVb(.  and  the  y  components  of  the  pressure  forces, 
such  that 
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F  =  pbc 


20, 


the  normalized  equation  for  the  transverse  direction  becomes 

(Fa  -  F^)  =  (pg  +  2Qg^)  sin  S 


(7b) 


Dividing  eq  (7b)  by  (7a)  and  assuming  that 

fPav  -  Pi)tic  tan  B  <<  (Pg  +  2Qg") 
the  two  equations  reduce  to 


(F= 


Fi) 


where 


tan 

»  -  .P.  * 

2Qs') 

(7c) 

Fa 

~  Pav*^c 

2Qca' 

(7d) 

2Qci' 

(7e) 

Fi 

=  Pibc  + 

b... 

The  supply  jet  deflection  is  thus  dependent  on  the  ratio  of  control-to- 
supply  momenta  and  pressures  forces.  An  indirect  dependence  on  the  off¬ 
set  D  and  the  wall  angle  a,  which  is  intuitively  suspected,  comes  in 
through  the  dependence  of  the  flows  and  interaction  region  pressures  on 
these  p ar ame te rs . 

2 . 4  Momentum  Equation  for  the  Attachment  Point 

The  momentum  flux  Jw  which  strikes  the  wall  at  an  angle  e  and  di¬ 
vides  at  the  returned  flow  streamline  nv  is  the  difference  between  the 
supply  momentum  J+  and  the  momentum  peeled  off  by  the  splitter  Jsp  (see 
fig.  2) .  Assuming  negligible  pressure  forces,  the  momentum  flux  balance 
equation  at  the  attachment  point  is 

Jw  ®  ~  Jdownstream  "  Jupstream 
or 


’w  (8) 

The  distances  ngp  and  n„  are  measured  from  the  jet  centerline  arc  to 
the  splitter  and  returned  flow  streamline  at  the  respective  jet  arc 
length.  Sgp  and  s^  (fig.  2  and  3).  The  first  integral  on  the  right 
in  eq  (8)  can  be  v.’rittcn  as 


ug^  -  uMn 

cos  0  = 

U>a„  -  1 

r  uMn 

J  -_oo 

S  =  Sgp_ 

_  lisp 

iiw 

fllw 
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\f  u^dn  t 

u'dn 

"sp 

L^"sp 

J  0 

and  since  the  momentum  between  any  two  lines  of  constant  similarity 
parameter  is  conserved 


rO 

fO 

u^n 

= 

uMn 

n 

s=s 

^  n 

sp 

w 

sp 

s=s 

sp 

Note  that  a(hgp/s)  is  constant  so  that  ngp  has  different  values  at 
s  =  s„  and  s  =  Sgp. 
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Using  the  velocity  profile,  eq  (3)  and  (4),  and  defining 
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and  noting  the  definition  of  from  eq  (6a)  and  (6b) ,  the  momentum 
flux  equation  degenerates  to 


-  T  +  T-  cos  0  -  -5-  T 


3  sp  sp  3 


T  _  -  I  +  2T  -  I-  T  5 


w  3  w 


The  root  of  the  cubic  equation  in  (eq  10)  applicable  to  the  attach¬ 
ment  point  is 


T  =2  cos 

V)  ( 


-  cos  ^  (X/2) 


X  =  1  +  I-  T  -  4  T  ’  + 
2  sp  2  sp 


fl  -  ^  T  +  4  T  M 
(  2  sp  2  sp  j 


The  Unknowns  in  3q  (11)  are  n„,  ngp,  s„,  Sgp ,  and  Except  for  the 

jet  virtual  origin  distance  the  unknowns  are  determined  through  the 

geometric  relationships  with  the  jet  position  in  the  amplifier. 


2 . 5  Upstream  Location  of  Supply  Apparent  Source,  syp 


Tile  previous  equations  involving  the  velocity  profile,  eq  (2),  con¬ 
tain  the  apparent  source  distance,  Sy^,.  To  determine  this  value,  the 
velocity  profile,  eq  (2) ,  is  equated  to  one-half  the  supply  nozzle  flow. 
Evaluating  at  s  =  0  results  in  the  expression  for  Sy^. 


and 
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(14b) 


2 . 6  Momentum  Flux  Equations  for  Jet  Curvature,  R 


An  arc  of  a  circle  is  assumed  for  the  attached  jet  centerline  curva¬ 
ture  (fig.  1  and  3) .  This  has  recently  been  experimentally  verified  for 
various  offsets  and  control  flows  (aspect  ratio  =  6,  turbulent  jet)  by 
Wada  and  Snimizu.^  The  jet  centerline  radius  of  curvature  is  the  ratio 
of  jet  momentum  flux  to  the  pressure  difference  across  it,  APj^  “  Pi  "  Pb 
so  that! 


R  = 


2Q  ^ 

or  ^p^  = 


(15) 


While  the  radius  of  curvature  of  the  jet  can  be  written  in  terms  of 
geometry,  both  Ap^,  and  Qs  in  eq  (15)  depend  on  the  bubble  pressure  p^ 
and  the  opposite-side  pressure  pj. 


Up  to  this  point,  the  difference  between  this  model  and  previous 
ones  (Epstein**  and  Lush,®  for  exeimple)  is  that  this  model  takes  into  ac¬ 
count  the  variation  in  flows  caused  by  changes  in  the  exit  pressures  and 
an  interaction  with  the  splitter.  Because  we  are  trying  to  describe  the 
two-sided  jet-attachment  device,  the  pressure  p^  on  the  opposite  side  of 
the  jet  is  not  assumed  to  be  zero  gage.  This  additional  unknown  re¬ 
quires  another  equation,  which  is  described  in  the  followd.ng  section. 


2 . 7  Vortex  Equation  for  the  Bubble  Pressure,  pb 


The  attachment  bubble  pressure  is  obtained  by  postulating  a  forced 
vortex  in  the  bubble.  For  simplicity  and  convenience,  it  is  assumed 
that  (fig.  4)  : 


(1)  The  piessure  distribution  is  given  by  the  differential 

equation  ^  =  2  * 

r  dr 

(where  m  is  the  angular  velocity  of  the  bubble  vortex) . 

(2)  The' driving  velocrty  at  the  edge  of  the  vortex  is  the 
entrainment  streamline  velocity,  Ug,  evaluated  at  Se  =  sw/2 
(fig.  4)  and  the  pressure  on  the  streamline  is  Pg . 

(3)  The  attachment  bubble  pressure  pi,  is  the  integrated 
average  of  the  vortex  pressure  distribution  resulting  in 


Pb  =  - 


+ 


P 


e 


(16) 
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2 . 8  Vortex  Driving  Velocity,  ue 


The  vortex  driving  velocity  Ue  is  obtained  from  the  velocity  pro¬ 
file  evaluated  at  the  downstream  distance  Se  and  cross-stream  distance 
He.  The  distance  s„,  and  hence  Se»  is  related  to  the  geometry,  and  the 
distance  n^  is  determined  by  integrating  the  velocity  profile,  eg  (2), 
between  n  =  0  and  n  =  n^  and  equating  the  flow  through  one-half  the 
supply  nozzle  which  results  in 


n 


(17) 


Substituting  eq  (17)  into  the  velocity  equation  (eq  2)  gives  the  vortex 
driving  velocity  as 
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Note  that  eq  (18)  is  not  accurate  near  the  supply  nozzle  exit  since 
the  entrainment  velocities  are  not  monotonically  decreasing  functions 
of  Sg  because  the  model  does  not  use  a  jet  potential  core  region.  Both 
the  turbulent  and  the  laminar  entrainment  streamline  velocities  increase 
to  a  maximum  and  then  decrease  for  increasing  values  of  Sg  in  the  model; 
whereas,  physically,  th;;  velocity  on  the  entrainment  streamline  is  high¬ 
est  near  the  nozzle  exit  and  decreases  as  the  downstream  distance 
increases . 

Rather  than  use  eq  (18)  in  the  region  where  Ug  is  inaccurate,  a 
cubic  distribution  correction  is  assumed  for  the  initial  downstream  dis¬ 
tances  to  insure  that  the  entrainment  velocity  monotonically  decreases 
from  the  supply  nozzle  velocity  Ug .  The  cubic  velocity  distribution  for 
Ug  is  matched  to  the  distribution  of  eq  (18)  at  the  distance  Sgm  associ¬ 
ated  with  the  maximum  velocity  Ugm.  With  the  derivatives  of  eq  (18) 
set  equal  to  zero,  the  maximum  Ue  and  its  location  are  determined. 
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The  cubic  function  for  s  <  s  is 

e  em 


Ug  =  niis^’  +  in2S^^  +  lOaSg  +  m4 


with  the  conditions 
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Solving  for  the  coefficients  (and  noting  that  u  =  Q  in  the  normalized 
form) 


mi 

'*s  em 

(21a) 

o  3 

®em 

ms 

=  -  2-  m,s 

2  ‘  em 

(21b) 

=  0 

(21c) 

=  Qs 

(21d) 

Equations  (18)  are  used  for  s  >  s  and  eq  (20)  and  (21)  for 
^  ,  e  em 

s  <  s 
e  em 


2.9  Geometric  Relations 


Tie  geometric  relations  presented  in  the  following  sections  are 
shown  in  figure  3.  The  attachment  wall  offset  distance  D  is  measured 
from  the  edge  of  the  supply  nozzle  to  the  beginning  of  the  attachment 
wall  plane,  and  the  attachment  wall  angle  a  is  measured  with  respect  to 
tile  x-axis  or  axis  of  symmetry.  The  vent  distance  dy,  the  jet  center- 
line  attachment  point  dw>  and  the  attachment  point  dap  are  all  measured 
from  the  exit  plane  of  the  supply  nozzle  along  the  attachment  wall. 

The  jet  radius  of  curvature  is  R.  The  centerline  arc  distance  s„  is 
measured  from  the  exit  plane  to  the  jet  centerline  attachment  point. 

The  deflection  angle  of  the  jet  is  f?  .  Both  the  jet  centerline  and 
attachment  streamline  (extended  to  the  wall)  intersect  the  wall  at  an 
angle  0  . 
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From  the  geometry,  the  jet  radi.us  of  curvature  is 


(D  4-  0.5)  cos  g 
cos  (a  +  65  -  cos  0 


the  jet  arc  length  is 

s„  =  R(c.  +  B  +  e) 

W 

and  the  jet  centerline  attachment  distance  is 

=  R[sin  (o(  +  6)  +  sin  9]  -  (D  +  0.5)sin  a 


(22) 


(23) 


(24) 


The  attachment  point  is  located  by  the  intersection  of  the  returned 
flow  streamline  with  the  wall.  Assuming  a  riglit-triangle  relationship 
with  the  jet  centerline  point  d^  and  attachment  angle  0  ,  the  attach¬ 
ment  point  is 

^ap  =  ‘^w  -  ®  ^25) 

(Vhen  the  jet  centerline  intersects  the  amplifier  axis  downstream  of  the 
splitter,  it  can  be  considered  as  a  switch,  since  more  than  half  of  the 
jet  flow  is  on  the  opposite  side.  The  intersection  of  the  jet  center- 
line  and  the  amplifier  axis  is 

Xj  =  2R  sin  B  (26) 


2 . 10  Attachment  Bubble  Volume  under  Jet  Centerline,  V 

The  volume  V  (per  unit  depth)  of  the  attachment  bxibble  is  bounded 
by  the  supply  nozzle  exit  plane,  the  attachment  wall  plane,  and  the 
jet  centerline  so  that 

V  =  0.5R^{a  +  6  +  6  -  0.5[sin  2 (a  +  6)  +  sin  29]} 

+  ER  sin(a  i  6)  -  0.5(D  +  0.5)  sin  a] (D  +  0.5)cos  a  (27) 

Because  o£  the  formulation  of  the  present  mJ^al,  there  is  the  pos¬ 
sibility  that  the  jet  centerline  may  intersect  the  opposite  wall. 

Rather  than  reformulate  all  of  the  equations  to  match  the  conditions 
imposed  by  the  solid  opposite  wall,  it  is  assumed  that  the  jet  curva¬ 
ture  remains  constant  even  though  it  would  pass  through  the  wall.  From 
figure  5,  the  jet  flow  corresponding  to  this  situation  is  assumed  to  be 
equivalent  to  the  jet  radius  of  curvature  R  deflected  Uirough  an  angle 
and  striking  the  opposite  wall  at  the  point  IN.  The  jet  centerline 
extends  through  the  wall  and  intersects  the  wall  again  at  point  OUT. 

In  reality,  however,  the  jet  flow  is  parallel  to  the  wall  along  the 
chord  length  IN-OUT.  At  OUT,  the  jet  flow  separates  from  the  wall 
and  continues  in  the  direction  of  the  jet  centerline.  Any  resultant 
wall  pressure  is  assumed  equivalent  to  the  pressure  required  to  change 
the  circular  path  of  the  jet  to  a  straight  path.  The  pressures  pi,  and 
Pl  are  affected  indirectly  only  by  the  different  volume  used  in  the 
calculations . 

The  attachment  bubble  is  reduced  by  the  area  of  the  segment  of 
circle  defined  by  the  points  IN  and  OUT.  Figure  5  shows  the  x-y  co¬ 
ordinate  system  with  the  origin  at  the  center  of  the  supply  nozzJe 
exit.  The  x-coordinates  of  the  points  IN  and  OUT  are  obtained  from 
the  solution  of  the  equations  for  the  opposite  wall 
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VOLUME  V 


Figure  5.  Jet  touching  opposite  wall 


y  »  X  tan  a  +  D  +  0.5 


and  the  jet  centerline  circle 


which  gives 


where 


=  (x  -  R  sin  g) *  +  (y  +  R  cos  8) 


xi  =  -  BW  -  /bWz  -  cw 
Xz  =  -  BW  +  »'BW2  -  CW 


BW  =  cos^  a  [(D  +  0.5  +  cos  6) tan  a  -  Rsin  6] 

CW  =  [ (D  +  0.5)cos  a] ^  +  2 (D  +  O.SjRcos  a  cos  8 


The  angle  i}i  subtending  the  arc  s^,  between  IN  and  OUT  is  given  by 

(j)  “  sin"^  jeln  8  -  +  sin~*  ^  -  sin  ^  (34) 

The  area  to  be  subtracted  from  eq  (27)  is 

=  0.5R^((|)  -  sin  (]))  (35) 

so  that  the  geometric  equation  for  the  attachment  bubble  volume  enclosed 
by  the  jet  centerline  is 


2.11  Lumped  Parameter  Line  Equations 


As  mentioned  previously,  the  supply,  control,  vent,  and  output 
channels  are  characterized  as  lumped-parameter  lines.  These  line 
equations  are  derived  by  integrating  Euler's  equation  along  a  stream- 
tube  assuming  one-dimensional,  incompressible  flow  with  no  body 
forces.  Along  some  streamtube  z,  the  normalized  integrated  form  of 
Euler's  equation  between  two  points  Zj  and  Zj  is 


Z2 


Zl 


9u 

at 


2  _ 


(37) 


Since  u(z,t)  =  Q(t)/A(z),  the  equation  can  be  written  as 


^  dz 

dt  A(z) 


+  U2^  -  Ul ^  =  Pi 


P2 


(38) 


The  integral  is  the  inductance  of  the  line,  but  for  convenience,  it  is 
multiplied  by  two,  so  that  the  inductance  parameter  L  is  defined  as 


L 


dz 

A(z) 


(39) 


The  inductance  parameter  depends  on  the  geometry  of  the  line  between 
any  point  Zi  chosen  as  the  input  (or  beginning)  of  the  line  and  the 
point  Zj,  which  terminates  inside  the  amplifier. 

The  remaining  terms  in  eq  (38)  Ccin  be  arranged  so  that  the  de¬ 
signer/experimenter  has  a  choice  of  using  either  the  static  or  the 
total  pressure  as  the  input  pressure  signal.  For  practical  applica¬ 
tion,  losses  must  be  accounted  for.  They  are  included  in  Euler's  equa¬ 
tion  by  using  the  discharge  coefficient,  c^  =  Qactual/Q ideal*  In  terms 
of  the  total  input  pressure,  p^.^  "*■  hj,  and  the  actual  flow  Q,  the 

line  equation  is 

3^  +  KQlQi  =  pti  -  P2  (40) 

where  the  resistance  parameter  K  is  defined  as 


K  H 


(43) 


where  Az  =  area  (per  unit  depth,  normalized)  at  station  zz.  Using 
q|q]  in  the  expression  allows  flow  in  a  line  to  be  in  either  direction. 
In  terms  of  static  pressure  input  pi ,  the  line  equation  is 


Cd 


+ 


(47) 


Of  course  when  (cdAz/Ai)^  <<  1  and  pi  =  Pti<  the  two  forms  of  the  line 
equations  become  equivalent.  The  discharge  coefficient  relationship 
is  of  the  form 


c^  =  (1  -  ai//^ie)  (1  -  az//^) 

where  ai  and  az  are  determined  from  the  solution  of  the  Karman-Pohl- 
hausen  momentum  integral  equations  for  a  particular  nozzle. 
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2.11.1  Supply  Equation 


The  normalized  supply  total  pressure  is  equal  to  one,  as  is  the 
normalized  supply  nozzle  width,  so  that  the  supply  resistance  param¬ 
eter  Kg  =  1/cjg^  .  The  supply  line  equation  is  then 


=ds 


+  K  Q  ^  =  1  - 
s  s 


(43) 


Note  that  no  reverse  flow  is  anticipated. 

Equation  (43)  couples  the  supply  flow  to  any  changes  in  the  interaction 
region  pressures. 

2.11.2  Control  Equations 

The  control  line  equations  are,  from  eq  (40) , 


dt 


do  • 
5t 
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^caQcalQcal 


KcAil^ci 


Pav 

Pi 


(44a) 


(44b) 


The  total  pressures  p^^,  Pd  are  the  input  forcing  functions  for  the 
entire  system  of  equations.  For  symmetrical  devices,  the  inductance 
parameter  Lc  ahd  the  discharge  coefficient  cjc  relationship  are  the  same 
for  both  control  lines.  The  resistance  for  the  controls  is  based  on  the 
distance  between  the  downstream  edge  of  the  control  and  the  edge  of  the 
jet.  The  edge  of  the  jet  is  determined  by  considering  the  jet  to  be  of 
constant  width  after  emerging  from  the  supply  nozzle  until  it  passes 
the  downstream  edge  of  the  controls.  At  the  downstream  edge  of  the  con¬ 
trol  (fig.  6,  distance  be),  the  amplifier  has  a  half-width  B.  The  ef¬ 
fective  opening  (tlirough  which  flow  may  pass)  of  the  attached-side  con¬ 
trol  is 

=  C-i  -  0.5)  +  Ay  (45a) 

ca  c 

where  Ay  is  the  distance  between  the  amplifier  axis  and  the  jet  cen¬ 
terline.  By  convention,  Ay  is  positive  when  i3  >  0 ,  The  opening  for 
the  opposite  control  is 


A^i  =  (B  -  0.5)  -  Ay^ 

The  resistance  of  the  controls  is  therefore 


(45b) 


K 


ca,ci 


1/A 


2 

ca  ,ci 


1/(B  -  0.5  ±  Ay^)^ 


(46) 


The  expression  for  Ay  .  is  obtained  from  the  equation  of  the  jet 
centerline  whose  radius  is 'R  and  whose  center  is  at  the  point  PJ  (fig.  6). 
The  circular  path  passes  through  the  point  PC(b^,Ay^)  so  tnat,  from  the 
equation  of  a  circle 


=  (b^  -  Rsin  B)*  +  (Ay^  +  Rcos  B)  ^ 

and  so 

Ay  =  -Rcos  B  +  (R^cos^  B  -  b  ^  +  2b  Rsin  B)*/^  (47) 

c  c  c 
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In  most  cases,  the  vent  ambient  pressures  are  equal,  so 
“  =  Pv  •  vent  line  equations  are,  from  cq  (40)  , 


<2Q., 


va  dt 


+  K  Q  Q 
va  va '  va ' 


=  P„ 


that 


dQ, 


VI 


vi  dt 


+  K  .Q  .  Q  . 

VI  VI '  vx ' 


=  Pv  '  Pi 


(50a) 

(50b) 


On  the  attached  side,  the  effective  vent  width  for  the  inductive  and 
resistive  parameters  is  based  on  the  distance  between  the  upstream 
edge  of  the  vent  and  the  relative  position  of  the  returned  flow  stream¬ 
line.  If  the  jet  attachment  point  is  upstream  of  the  vent,  the  vent 
flow  into  the  bubble  is  zero.  As  tile  attachment  point  sweeps  past  the 
vent,  the  effective  opening  increases  from  zero  to  the  phy.sical  width 
of  the  vent  channel.  The  attached-side  vent  inductive  parameter  Lva 
is  variable  to  compensate  in  some  measure  for  the  reduction  in  vent 
area  caused  by  the  jet  flow  that  is  directed  into  the  vent  as  the  jet 
centerline  is  passing  the  vent.  The  minimum  resistance  K^min  based  on 
the  mascimum  vent  opening  is  then 


K 


vmin 


1/W 


(51) 


On  the  opposite  side,  the  vent  is  never  covered  by  the  jet  (at 
least  for  the  anticipated  geometries  considered  here);  therefore, 

Kvi  ~  Kvmln  =  constant.  On  the  attached  side,  however,  resistance 
values  greater  than  Kvmln  occur  depending  on  the  effective  vent  opening 
Ava  (fig*  6) •  Tnis  distance  is  measured  along  the  jet  radial  line  pass¬ 
ing  through  the  upstream  corner  of  the  vent  and  extending  to  the  attach¬ 
ment  streamline.  The  arc  length  Sy  of  the  jet  centerline  to  the  radial 
line  passing  through  the  vent  corner  is 


=  R(6  +  C) 

where 


(52) 


tan" ' 


d^cos  a  -  R  Cos  6 
_R  cos  6  -  (D  +  b .5)  -  s^sin  a 


(53) 


At  the  arc  distance  Sv,  the  distance  from  the  jet  centerline  to  the 
attachment  streamline  is  nv.  The  expression  for  nv  is  obtained  by  using 
the  similarity  properties  for  the  jet  solution.  Since  the  similarity 
parameter  is  constant  for  a  particular  streamline,  in  this  case,  the  re¬ 
turned  flow  streamline,  its  value  is  obtained  at  the  jet  centerJ.ine 
attachment  point,  s  =  s„.  The  value  of  the  similarity  parameters  for 
tile  at taclimerit  streamline  are 

on 

n  =  —  TURBULENT  (54a) 


5  =  — TT—  laminar  (54b) 

\  w 

where  the  n  's  are  gWen  by  eq  (12).  The  distance  n  is  obtained  by 
combining  eq  (52)  th- ough  (54)  yielding 
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The  effective  restriction  for  vent  flow  is  therefore 


(55b) 


R 


'v  ■■  ^v  sin  c 

and  the  resistance  parameter  for  the  vent  is 

K  =  1/A  ^ 

V  v 


d  cos  a  -  R  sin 

V 


(56) 


(57) 


when  A  <  W  .  The  area  A  is  then  also  used  in  the  inductive  param¬ 
eter  wheX  the  jet  is  sweeping  by  the  vent. 

2.11,4  Output  Equations 

The  equations  for  the  output  lines  are,  from  aq  (40), 
dQ 

L  +  K  Q  IQ  I  =  (p-+Pj)-P 

oa  dt  oa  oa '  oa '  *^1  ^da  ^oa 


(58a) 


dQ. 


^oi  3r 


+  K  .Q  .  Q  . 

01  OL ’  01 


i!  =  (Pi  +  Pai)  -  p. 


(58b) 


The  pressures  p^^  and  p^^  are  the  static  pressures  at  the  external 
measuring  point  of  the  output  lines.  The  total  pressures  (p^  +  Pda) 
and  (pi  +  Pdi)  provide  the  internal  driving  force  to  the  output  lines. 
Since  the  active  output  line  is  for  the  most  part  separated  from  the 
bubble  by  the  jet,  the  static  pressure  at  the  entrance  to  it,  as  well 
as  to  the  inactive  output  line,  is  p^ . 

The  dynamic  pressures  p^jg  and  p^j^  are  determined  from  the  momentum 
flux  impinging  on  the  output  receivers  of  area  Woa  =  WqI  =  Wq.  The 
momentum  flux  impinging  on  the  active  output  line  is  that  which  is  flow¬ 
ing  dowaatream  from  the  attachment  point.  Averaging  the  momentum  over 
the  entrance  area  of  the  output  line  yields  the  dynamic  pressure  such 
that : 


■da 


downstream 


(59) 


where  Jdownstream  is  given  by  eq  (8) .  Equation  (8)  also  includes  the 
momentum  peeled  off  by  the  splitter.  Averaging  this  momentum  over  the 
inactive  output  line  opening  gives  the  dynamic  pressure  pj  impinging 
on  that  output  as 

p^.  =  i  splitter  (bO) 

Ml  2  W 


The  dynarr.ic  pressures  can  be  expressed  in  terms  of  the  attachment  point 
parameters  Tsp  and  T'v  (eg  6,  9,  10,  and  11). 


Pda 

^di 


0.75Q  Mt  h  T  -  (T  ^  +  T  ^)/3]/W 
sp  w  sp  w  '  '  o 

0.75Q  '  [2/3  "  T  +  T  V31/W 
s  sp  sp  '  o 


(61a) 

(61b) 
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Tlie  resistance  parameters  of  *-he  outputs  are  constant 

at  Koa  =  1^0 i  "  inductive  parameters  of  all  lines  are 

a  function  of  their  shape  and  are  obtained  by  using  eq  (39), 

Various  expressions  have  been  developed  nerein  to  describe  the  bi¬ 
stable  amplifier  model.  Tne  line  equations  couple  the  conditions  ex¬ 
ternal  to  the  amplifier  to  the  internal  conditions;  however,  the  coup¬ 
ling  of  the  equations  in  the  internal  region  of  the  amplifier  depends 
on  the  pressure  level  term  Pg  (eq  16).  A  constant  control  volume  around 
the  entire  amplifier  is  used  to  calculate  the  pressure  p^,  which  is  con¬ 
tinuously  adjusted  so  that  the  net  flow  is  zero.  Thus,  since  the  rela¬ 
tionships  for  all  the  into  the  unit  are  known  as  functions  of  the 

respective  pressure  differences,  there  is  always  one  value  of  p^  (and 
hence  of  all  the  other  pressures)  that  will  satisfy  continuity  at  any 
given  time.  By  calculating  all  the  line  flows  (as  determined  by  the 
respective  pressures)  and  summing  them,  a  non-zero  resu].t  may  occur. 

If  it  does,  tnen  p^  is  adjusted  so  that  the  net  flow  is  zero. 


3  .  .IBTHOD  OF  SOLUTION 

The  geometric  and  flow  equations  (1  through  61)  are  solved,  simul¬ 
taneously  where  necessary,  on  a  digital  computer.  The  programming 
language  is  the  Digital  Simulation  Language  DSt/go,'^  a  sub-language 
based  on  FORTRAH  IV  with  a  large  built-in  library  of  internal  functions 
and  subroutines  which  allow  simple  commands  to  be  used  when  solving 
coraple,:.  differential  and  implicit  equations.  Of  the  several  integration 
routines  available  in  DSL/90,  the  technique  chosen  was  Milne,  a  fifth- 
order,  predictor-corrector  routine  with  variable  time-step  size. 

Equations  (4),  (43),  (44),  (50),  and  (58)  are  integrated  in  time. 
Arbitrary,  but  reasonable,  values  are  assinged  as  initial  conditions 
for  the  integrals.  The  program  then  takes  the  assigned  values  of 
the  program,  geometric  and  fluid  constants  and  parcumeters ,  and  in  con- 
juriction  with  the  governing  set  of  equations,  relaxes  the  problem 
variables  to  the  actual  steady-state  values.  Tliese  values  then  consti¬ 
tute  tlie  initial  conditions  for  the  transient  response  of  the  bistable 
amplifier  to  an  arbitrary  input  pressure  signal  which  is  impressed  after 
the  initial  relaxation  time. 

The  program  listings  are  presented  in  appendix  A.  The  program  names 
and  variables  are  defined  tliere  as  are  the  specific  instructions  for 
program  start-up.  In  the  calculation  of  the  incremental  change  in  pres¬ 
sure  level,  an  iterative  procedure  at  each  time  step  is  specified  during 
tlie  simultaneous  solution  calculation  of  the  entire  set  of  equations. 
Subsequent  iterations  are  based  on  the  new  pressure  level.  This  itera¬ 
tion  at  each  time  step  used  excessive  machine  time  and  so  was  discarded 
in  favor  of  the  following  alternate  method.  The  sum  of  the  flows  cal¬ 
culated  from  the  various  equations  was  generally  non-zero.  Instead  of 
Iterating  to  find  a  true  value  of  p^  at  a  particular  time  step,  a  change 
was  calculated  simply  as  a  function  i  C  the  net  flow.  The  error  incurred 
was  less  than  2  percent,  overall. 


1  2 


Syn,  W.M.  and  Wyraan,  D.G.,  "Digital  Simulation  Language  User's  Guide,"  IBM  Systems 
Development  Division,  San  Jose,  California,  TR  02.355,  01  July  1965, 
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4. 


liESULTS 


Three  levels  of  complexity  are  presented.  First,  the  theory  in  its 
simplest  form  is  used  to  predict  the  switching  dynamics  of  fluid  ampli¬ 
fiers  for  the  case  where  splitter  and  output  effects  are  negligible. 

An  intermediate  formulation  is  used  to  predict  the  switch  time  (time 
to  reach  the  splitter)  of  the  jet  when  the  splitter  and  output  effects 
are  not  .legligible.  In  this  case,  however,  the  jet  remains  curved 
towards  the  initial  attachment  wall.  Finally,  the  most  complex  model 
is  used  to  predict  the  response  of  an  entire  fluid  amplifier,  in  which 
the  jet  is  allowed  to  completely  reattach  to  the  opposite  side  by  allow¬ 
ing  a  reversal  in  jet  curvature.  The  physical  constants  for  air  were 
used.  A  constant  value  of  10  was  used  throughout  for  the  spread  param¬ 
eter,  a  . 

4 . 1  Comparison  of  Turbulent  Attached  Jet  Results  with  Published  Data 

Before  initiating  the  Harry  Diamond  Laboratories  (HDL)  experimental 
effort,  the  simple  analytical  results  were  f..rst  compared  to  existing 
published  data.  These  results  may  also  be  found  in  Goto  and  Drzewiecki.’ 

4.1.1  Steady-State  Attached  Turbulent  Jet,  Negligible  Splitter 
and  Output  Effects  (Simplest  Model) 

The  steady-state  attached  jet  solution  considers  a  two-wall  geometry 
with  closed  control  ports.  With  no  control  flow,  there  is  no  momentum 
deflection  of  the  jet  and  the  effect  of  the  opposite  wall  is  minimized 
although  not  removed.'  ^  The  essential  difference  between  this  case  and 
the  single-wall,  steady-state  model  is  that  for  a  given  supply  pressure, 
the  two-wall  model  has  more  supply  flow  because  of  the  negative  pressure 
in  the  interaction  region.  Figure  7  shows  the  agreement  between  the  pre¬ 
dicted  results  and  the  experimental  data  of  Kimura  and  Mitsuoka'  for  the 
attachment  distance  dap  versus  offset  D  for  a  wall  angle  of  0.262  rad 
(15  deg),  with  aspect  ratio  (AR)  of  about  3.3. 

Lush®  (fig.  VIII,  p.  36)  gives  steady-state  experimental  data  for 
turbulent  jet  attachment  distance  versus  offset  for  open  control  ports. 
The  analytical  results  for  Lush's  geometry  were  obtained  by  opening  both 
controls  to  zero  (gage)  input  pressure  at  time  zero  and  allowing  the  jet 
to  reach  a  new  equilib  ium  position.  As  shown  in  figure  7,  the  agree¬ 
ment  between  the  present  theory  and  Lush's  experimental  data  is  good  for 
the  given  range,  0  f  D  <  1.0,  for  a  wall  angle  of  0.262  rad  (15  deg)  and 
AR  =  1.0. 


4,1.2  Switch  Time  of  the  Turbulent  Attached  Jet  with  Negligible 
Splitter  and  Output  Effects 

The  experimental  switch  time  data  of  Johnston''*  and  Lush®  are  com¬ 
pared  here  with  theoretical  results.  Johnston's  data  was  obtained  on 
two  test  models  that  had  no  splitter  or  side  wall  vents.  One  unit  had 
a  length  of  nine  nozzle  widths,  the  other  17.  The  dimensions  are  given 
by  Johnston'"  and  also  by  Goto  and  Drzewiecki®  and  are  shown  in  figures 
8  and  9 . 


"Perry,  C.C,,  "Two-Dimensional  Jet  Attachment,"  Advances  in  Pluidico,  Proc,  ASMP 
Fluidics  Symposium,  Chicago,  May  1967. 

'"Johnston,  R.P.,  "Dynamic  Studies  of  Turbulent  Keattachraent  Fluid  \mpllfiers," 
Master's  Thesis,  University  of  Pittsburgh,  School  of  Engineering  and  Mines,  1963, 
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Figure  7-  Turbulent  attachment  distunce  versus  offset 


The  time  it  takes  for  the  jet  centerline  to  intersect  the  unit 
axis  at  9.0  or  17.0  nozzle  v/idths,  the  ends  of  the  units,  is  considered 
to  be  the  switching  time.  Switching  times  were  calculated  using  a 
signal  that  reached  the  selected  pressure  after  a  200-  us  ramp  to  dup¬ 
licate  Johnston's  input  signal.  Figures  8  and  9  show  the  experimental 
data  and  the  analytical  results  of  switching  time  versus  the  various 
final  input  pressure  amplitudes.  For  each  case,  both  the  analytical 
curve  for  the  attachment  point  reaching  the  end  of  the  wall  (separation) 
and  the  jet  centerline  intersecting  the  axis  at  the  end  of  the  model 
are  shown. 

The  agieement  between  theoretical  switching  time  and  experimental 
data  is  better  for  the  short-wall  model  (fig.  8)  than  fer  the  long-wall 
model  (fig.  9).  In  the  latter  case,  at  the  lower  values  of  input 
pressure,  the  jet  radius  of  cu  vature  is  quite  large,  so  that  when  the 
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re  8.  Turbulent  jet  switching  time  versus  control  pressure,  0.2  ms 
ramp  input,  Johnston's  short  wall  data 


Figure  9.  Turbulent  Jet  switching  time  versus  control  pressure,  0.2  ms 


ramp  input,  Johnston's  long  wall  da^a 


centerline  intersects  the  axis  at  17  nozzle  widths,  a  major  portion  of 
the  jet  and  the  vortex  are  outside  the  model  geometry.  The  analytical 
model  no  longer  applies  when  the  jet  is  outside  the  geometry,  since  one 
no  longer  expects  tne  jet  to  curve.  The  effect  of  such  a  lack  of  curva¬ 
ture  would  be  an  increase  in  vent  flow  eind,  therefore,  a  decrease  in 
switch  time.  In  general,  the  agreement  for  both  geometries  is  fairly 
good,  since  the  line  dimensions  and  signal  characteristics  were  esti¬ 
mated  from  photographs  presented  by  Johnston.'*' 

Lush's  test  model  had  a  splitter  at  20  nozzie  widths  downstream 
(to  minimize  splitter  effects)  and  a  0.262-radian  (15-deg)  wall 
angle.  The  supply  and  control  nozzles  were  2-54  cm  (1.0  in.)  with  an 
aspect  ratio  of  one.  Offset  was  0.482,  dv  was  13.04,  and  Ay  was  2.2. 

The  other  dimensions  were  estimated  and,  in  nozzle  widths,  are: 
the  supply  chamber  length,  5;  width,  3;  control  chamber  length,  3; 
width,  1;  and  vent  length,  3.  Details  of  his  geometry  are  given 
by  Lush,®  In  his  model,  the  controls  were  initially  open  to  atmosphere. 
His  10-  to  20-ms  pressure  rise-time  input  signal  (Lush,®  p.  76)  is 
represented  in  this  model  by  a  20-ms  ramp  to  a  final  pressure  amplitude. 
The  experimental  switch  time  was  measured  from  the  beginning  of  the 
pressure-flow  rise  in  the  control  chamber  to  the  time  the  jet  reached 
the  end  of  the  opposite  wall.  Using  the  initial  simple  theory,  the 
analytical  switch  time  was  started  at  the  same  point  but  ended  when  the 
jet  centerline  intersected  the  amplifier  axis  at  xgp  =  20  (the  splitter 
distance) .  Figure  10  shows  the  experimental  and  analytical  results 
of  switching  time  versus  final  control-pressure  amplitude.  As  is  shown, 
the  theory  agrees  well  with  the  data. 

No  comparison  has  been  made  with  Muller's  data,'®  since  there  was 
not  sufficient  information  about  his  splitter  shape  and  location,  out¬ 
put  line  dimensions,  and  input  measuring  locations.  His  device  was  a 
symmetrical,  low-offset  device,  so  it  was  unfortunate  that  comparisons 
could  not  be  made, 

4 . 2  HDL  Experimental  Apparatus  and  Methods 

The  HDL  experimental  program  was  initiated  after  obtaining  good 
agreement  between  the  initial  simple  analysis  and  existing  published 
data.  The  purpose  of  the  experimental  program  was  to  compare  the 
analyses  of  the  simple  model,  the  intermediate  model,  and  the  most 
complex  model  to  data  of  amplifier  characteristics  such  as  steady-state 
recovery,  switching  characteristics,  and  actual  amplifier  response 
time . 

4,2.1  Instrumentation 

Figure  11  is  a  diagram  and  figure  12  a  photograph  of  the  HDL  test 
set-up.  A  not-film  anemometer  probe  was  flush-mounted  in  the  test 
model  cover  plate  and  located  at  a  given  point  on  the  geometric  axis 
of  the  device.  As  the  jet  switched,  the  maximum,  or  centerline  velocity 
passed  that  point,  and  the  anemometer  registered  a  maximum  voltage. 

The  response  of  the  constant-temperature ,  hot-film  system  was  greater 
than  20  kHz;  however,  the  output  signal  was  electronically  filtered  to 
pass  only  5  kHz.  The  control  signal  was  initiated  by  an  electromechani¬ 
cal  solenoid  valve  upstream  from  the  control  port  connection  to  the 
test  model.  The  control  flow  entered  tlie  device  perpendicular  to  the 


'®Huller,  H.R.,  "Wall.  Re  attachment  Device  with  Pulped  Control  Flow,"  Proc.  2nd  Fluid 
Amplification  Symposium,  HDL,  Vol.  I,  May  1964. 
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figure  10.  Turbulent  jet  switching  time  versus  control  pressure,  2C 


ramp  input,  Lush's  data 
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Figure  11.  HDL  test  setup. 


cover  plate  and  impinged  directly  on  a  strain-gage  type  pressure  trans¬ 
ducer  mounted  in  the  bottom  of  the  control  channel.  T.he  transducer 
measured  the  total  pressure  of  the  irapinging  flow.  (This  transducer 
measurement  of  total  pressure  was  verified  by  actual  total  pressure 
probe  measurements.)  The  response  of  the  pressure  measuring  system 
was  less  than  0.05  ms  and  was  neglected,  since  the  rise  time  of  the 
input  signal  was  around  1.5  ms.  Output  total  pressures  were  measured 
with  a  pitot  probe  which  had  a  3-inm  diameter  and  was  2  cm  long.  It  was 
fitted  directly  into  a  strain-gage  type  pressure  transducer.  The  re¬ 
sponse  was  estimated  at  over  4  kHz.  The  control  signal  and  the  hot- 
film  output  were  monitors  1  simultaneously  on  a  dual-trace  oscilloscope 
to  determine  the  time  hi  tory  of  the  control  signal  and  the  time  to  move 
the  jet  centerline  past  .he  hot-film.  Steady  state  control  flows  into 
the  amplifier  configuration  were  measured  with  laminar-flow  meters, 

4,2.2  Experimental  Models 

Table  I  lists  the  dimensions  of  the  three  large  configurations 
tested.  Except  for  the  offsets  D,  control  line  areas  Wc ,  ,ind  the 
depth.s ,  the  models  had  the  Scime  plan -view  geometry.  The  splitter, 
when  used,  was  a  plexiglas  wedge  with  an  included  angle  of  0.418  rad 
(24  deg)  that  could  be  located  at  any  desired  position  in  the  configura¬ 
tion.  The  location  cf  the  splitter  was  chosen  as  10  nozzle  widths  down¬ 
stream  of  the  power  nozzle,  since  splitters  in  most  actual  fluidic  de¬ 
vices  are  located  at  about  that  distance. 
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LINE  SIZES 

Supply : 

Vent  : 

length 

=  30 

rum 

length 

=  20 

mm 

width , 

w  =. 

10  mm 

width , 

W  = 

4  I 

Control : 

Output : 

length 

-  20 

mm 

length 

=  56 

mm 

width. 

Wr  = 

4  mnr  (20  nuii  for 

Model  3) 

width , 

W  *  = 

7  1 

The  pattern  for  experimental  Models  1  and  2  was  programmed  and  then 
drawn  by  a  digital  computer  in  conjunction  with  a  Cal-Comp  plotter 
using  the  AUTOSKKM  program  shown  in  appendix  B.  The  plotter  also  pro¬ 
vided  a  twice  size  {2:1)  cut-and-strip  negative  of  the  device.  Using  a 
photoetching  process  (Dycril) ,  a  template  in  plastic  was  obtained  so 
that  the  models  could  be  manufactured  on  a  pemtograph  milling  machine. 
Model  1  was  milled  in  aluminum,  and  Model  2  in  phenolic  plastic.  The 
control  line  cross  sectional  areas  of  Model  2  were  later  enlarged,  re¬ 
sulting  in  Model  3, 

4.2.3  Discharge  Coefficients,  cj 

Instead  of  solving  for  the  discharge  coefficient  c^  (sec  2.11),  the 
discharge  characteristics  of  the  supply  and  control  nozzles  in  this 
study  v/ere  determined  experimentally.  As  previously  noted,  the  nozzles 
for  all  models  were  the  same  width.  The  supply  nozzles  had  a  contrac¬ 
tion  ratio  (cross  sectional  area  of  line  to  orifice  area)  for  Models  1 
and  2  of  six  and  a  contraction  length  ratio  (length  over  which  contrac¬ 
tion  occurs  to  orifice  area)  of  five.  The  control  nozzles  for  Models 
1  and  2  had  a  contraction  ratio  of  two  and  a  contraction  length  ratio 
of  0,875.  Modifying  the  Model  2  control  line  resulted  in  a  contraction 
ratio  of  10  (Model  3) .  Figure  13  shows  the  supply  and  control  geometry 
for  Model  2,  which  is  representative  of  all  tJie  models  (except  for  I’;. 
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Figure  14  shows  the  discharge  coefficients  cj  as  a  function  of 
Reynolds  number,  Q'/h'v'  (based  on  nozzle  width  and  actual  flow),  for 
the  test  models.  The  supply  cja's  were  the  same  for  all  models  within 
the  experimental  error.  Two  discharge  curves  were  obtained  for  the  con¬ 
trol  nozzles.  The  upper  curve  for  the  control  discharge  coefficient  cjc 
is  the  actual  cdc  curve  for  the  nozzles.  The  lower  curve  is  actually 
a  line-loss  coefficient,  because  the  point  of  measurement  upstream  from 
the  nozzle  included  an  elbow  fitting  as  well  as  the  control  nozzle. 

The  curve  is  presented  because  some  switch  time  experiments  (Model  1 
without  the  splitter)  had  used  that  measurement  point  to  monitor  the 
control  input  signal. 

For  simplicity  in  the  analysis,  the  supply  discharge  coefficient 
was  assumed  constant  because  the  variation  of  supply  flow  during  a 
switching  calculation  was  small  compared  to  the  variation  of  cds  with 
Re.  The  turbulent  switching  experiments  were  run  at  Reynolds  numbers 
of  Re4.  =  15,000  to  18,000,  and  the  laminar  at  Re+  -  1000.  Supply  dis¬ 
charge  coefficients  of  Cds  =  0.85  and  0.80,  respectively,  were  used. 
Control  coefficients  for  the  analysis  were  based  on  a  linear  approxi¬ 
mation  to  the  data,  because  the  control  flow  and  hence,  control  Rey¬ 
nolds  number  varied  significantly  during  the  time-dependent  calculations. 

4 . 3  Comparison  of  Turbulent  Attached  Jet  Results  with  HDL  Data 

4.3.1  Switching  Time  of  the  Turbulent  Attached  Jet  with  Negligible 
Splitter  and  Output  Effects 

.Model  1  was  used  to  obtain  the  time  response  of  a  turbulent,  at¬ 
tached  jet  without  splitter  effects.  The  supply  pressure  was  fixed  at 
7,5  kPa  (1.09  psig) .  The  input-pressure  signal  was  a  ramp  to  a  pre¬ 
selected  level  and  was  applied  to  the  attached-side  control.  The  ramp 
duration  from  the  initial  pressure,  Pgvf  final  input  value  p^^ 

varied  from  between  1  and  2  ms .  In  the  analytical  computations,  the 
input  ramp  rise-time  was  fixed  at  1.5  ms. 

34 


The  hot-film  probe  was  located  at  10  nozzle  widths  downstream.  The 
time  at  which  the  probe  registered  rhe  maximum  signal  was  used  as  the 
experimental  switch  time.  The  time  at  which  the  jet  centerline  inter¬ 
sected  the  amplifier  axis  at  10  nozzle  widths  was  used  as  the  analytical 
switch  time.  Control  pressure  versus  switch  time  data  was  collected  for 
two  conditions  of  inactive  control  port  loading:  blocked,  and  open  to 
ambient  pressure.  Figure  15  shows  the  data  with  the  inactive  control 
blocked  (Qd  =  0) ,  and  the  agreement  is  excellent  over  the  entire  range 
of  input  control  pressures.  Figure  16  shows  the  data  where  the  inactive 
control  was  open  to  ambisnt  pressure  (Qci  /■  0).  The  agreement  is  fairly 
good  over  the  range,  but  the  theory  slightly  underestimates  the  switch 
time  at  the  higher  input  pressures.  As  expected,  both  the  theoretical 
and  experimental  switch  times  are  slower  if  the  inactive  control  is  open 
rather  than  blocked,  because  there  is  a  switch-retarding  flow  coming  in 
through  the  inactive  control,  making  the  opposite  control  pressure  higher. 

4.3.2  Switch  Time  of  the  Turbulent  Attached  Jet  with  Splitter  and 
Output  Effects  (Intermediate  Model) 

Experimental  data  were  obtained  on  two  flueric  amplifiers.  Models  1 
and  2.  Both  had  sharp  splitters  with  their  leading  edges  located  at  10 
nozzle  widths  downstream  of  the  supply  exit  plane  where  b4.  =  2  miti.  The 
supply  pressures  were  10  kPa  and  u.6  kPa,  respectively.  Again,  input- 
pressure  signals  consisting  of  7.5-ms  ramps  to  a  final  pre-selected 
level  were  used.  Figures  17-20  show  the  comparison  between  theory  and 
data.  With  exception  of  the  case  in  figure  19,  the  agreement  is  good. 

No  immediate  explanation  is  offered  for  the  discrepancy  at  low  values 
of  input  pressure  shown  in  figure  19. 

4.3.3  Response  Time  of  a  Flueric  Wall-Attachment  Amplifier 
(Intermediate  Model) 


The  total  pressures  of  the  outputs  during  switching  were  measured 
with  a  pitot  probe  located  at  the  exit  plane  of  the  output  line.  Meas¬ 
urement  of  both  the  active  and  inactive  outputs  ddring  switching  were 
made,  but  not  simultaneously.  Model  3  was  used  in  this  test,  with  the 
supply  pressure  at  7.5  kPa.  The  inactive  control  was  blocked  and  the 
input  signal  final  amplitude  was  0.25  of  the  supply  pressure,  after  a 
1.5 -ms  duration  ramp. 

Figure  21  shows  the  comparison  between  experiment  and  theory.  The 
agreement  is  good  in  the  region  up  to  and  just  past  the  point  where  the 
jet  centerline  passes  the  splitter  (txj=io);  past  that  point,  the  data 
and  theory  diverge  because  the  intermediate  theory  does  not  allow  for  a 
reversal  of  jet  curvature.  Just  after  the  jet  centerline  passes  the 
point  of  the  splitter,  the  jet  actually  attaches  to  the  opposite  side, 
whereas,  the  intermediate  model  centerline  attachment  point  continues 
to  sweep  downstream  on  the  original  attachment  side. 

In  general,  if  one  considers  the  theoretical  amplifier  response  to 
be  the  time  of  intersection  of  the  active  and  inactive  pressure  traces, 
then  agreement  with  experiment  is  very  good  when  compared  to  the  rise 
time  of  the  inactive  output  signal.  Figure  21  shows  that  the  inactive 
output  reaches  90  percent  of  its  final  value  between  5  and  5.5  ms,  and 
the  theoretical  output  curves  intersect  at  about  5.5  ms. 

4.3.4  Steady  State  Characteristics  of  Flueric  Wall  Attachment 
Ampiltiers  Unrermediate  Model) 

As  a  further  test  of  the  design  utility  of  the  analytical  model, 
several  theoretical  and  experimental  steady  state  values  were  compared. 

35 


normalized  control  pressure  .  p„  ,  ,  normalized  control^  pressure 


Figure  15.  Control  pressure  amplitude  versus  switch  time — no  splitter, 
inactive  control  blocked,  D  =  0,9. 


Figure  16.  Control  pressure  amplitude  versus  switch  time — no  splitter, 
Inactive  control  blocked,  D  =*  0.9. 
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NORMALIZED  CONTROL  PRESSURE 


MODEL  2  WITH  SPLITTER  INACTIVE  CONTROL  fitOCKED  p+  *  T5MP 


Figure  21.  Transient  output  total  pressure  versus  time 


The  control  characteristics  Qj,  versus  P,.  were  obtained  from  each  side 
of  the  amplifier  lor  both  open  and  blocked  inactive  control  ports.  Al¬ 
though  the  characteristics  for  the  two  controls  of  a  model  were  slightly 
different,  they  were  averaged  to  maintain  clarity  (fig.  22  and  23) .  The 
points  of  comparison  with  theory  are  the  pressure  at  zero  flow  and  the 
pressure  and  flow  at  switching. 
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The  value  of  the  initial  condition  of  the  active  control  exit  pres¬ 
sure  (Pav  at  t  =  0)  of  the  present  tests  corresponds  to  the  experimental 
pressure  at  zero  flow.  The  value  of  minimum  pressure  to  switch  corres¬ 
ponds  to  the  experimental  steady-state  switch  pressure.  Theoretical 
values  of  minimum  pressure  to  switch  were  not  determined  explicitly 
since  the  control  pressure  was  pre-selected  in  the  computer  program. 

For  this  reason,  a  range  of  values  for  the  theoretical  minimum  switch 
pressure  is  used  for  the  comparison  with  experiments.  Tlie  theoretical 
values  shown  in  figures  22  and  23  show  that  the  agreement  is  good. 

The  experimental,  dynamic-pressure  transfer  characteristics  for  the 
active  output,  Po (=  Pi  +  Pda^  versus  p^a  shown  in  figures  24  and  25. 

The  outputs  are  completely  open  so  that  the  output  pressure  is  the  total 
pressure  measured  with  a  pitot  probe  at  the  exit  plane  of  the  output 
lines.  Because  of  some  slight  asymmetry  encountered,  transfer  curves 
for  both  outputs  of  Model  1  are  given,  whereas  a  single  curve  is  repre¬ 
sentative  of  measurements  on  either  side  of  Model  2.  The  points  of  com¬ 
parison  with  theory  are  the  total  pressures  of  the  active  output  when 
the  active  control  is  blocked,  open  to  ambient,  and  at  the  switching 
pressure  level.  The  respective  theoretical  quantities  are  determined 
from  the  values  of  the  flow  out  of  the  active  output  converted  to  a  dy¬ 
namic  pressure  when:  the  control  input  pressure  level  equals  the  control 
exit  pressiire  at  time  zero  (Pav  t  =  0);  the  input  pressure  signal 
level  is  zero  gage;  and  at  the  minimum  switching  pressure  level.  The 
comparison  is  made  in  figures  24  and  25  with  the  agreement  being  good  in 
most  respects. 

The  comparison  of  experimental  and  analytical  steady  state  values 
associated  with  the  characteristics  curves  presented  above  are  summar¬ 
ized  in  table  II,  which  lists  their  numerical  values. 


4.3.5  Miscellaneous  Results 

Since  the  analytical  and  experimental  results  obtained  agree  favor¬ 
ably,  unverified  analytical  results  can  be  presented  with  some  measure 
of  confidence.  Results  for  any  variable  can  be  obtained  from  the  compu¬ 
ter  program,  but  the  attachment  point,  active-control  flov/,  and  the  in¬ 
active-output  flow  are  chosen  as  illustrative  of  the  events  during  a 
switch.  Figure  2G  shows  these  variables  as  a  function  of  the  dimension¬ 
less  time  (t  =  t'/t+)  real  time  t'.  The  results  correspoi’d  to  the 

conditions  shown  in  figure  19  for  Model  2  with  a  splitter,  p|  =  7.5  kPa, 
inactive  control  blocked,  and  p^g  0.1. 

The  attachment  point  distance  dgp  and  the  control  flov?  Qca  increase 
rapidly,  until  the  final  control  pressure  p^g  =  0.1  is  reached  at  t'  = 
1.5  ms.  Note  that  Qca  has  a  small  overshoot.  Although  the  control  pres¬ 
sure  is  now  constant,  the  bubble  pressure  is  changing  and  the  attachment 
point  distance  and  control  flow  once  again  increase  until  the  attachm.ent 
point  reaches  the  vent  (d^  =  9) .  Uncovering  the  vent  changes  the  con¬ 
trol  flow  into  the  bubble  so  that  it  fluctuates  while  the  attachment 
point  dwells  on  the  edge  of  the  vent.  The  attachment  point  distance 
then  increases  again,  jet  switching  starts,  and  the  splitter  begins  to 
intercept  flow.  This  is  shown  by  the  increase  of  the  inactive  output 
flow  Qoi  when  switching  begins.  The  motion  of  the  attachment  point 
si.ows  down  on  "seeing"  that  output's  inductance.  Muller’ ^  shows  control 
flow  overshoot  and  slight  oscillation  in  much  the  same  manner  as  does 

the  pre.sent  analysis. 
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Figure  24.  Transfer  characteristics.  Model  1. 
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Figure  25.  Transfer  characteristics.  Model  2. 
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Figure  26.  Analytical  turbulent  output  variables  during  switching,  Model  2. 


Table  II.  Steady  State  Characteristics, 


Steady  State  Dy.iamic  Pressure  Recovery  (pj  4- 
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4 . 4  Comparison  of  Laminar  Attached  Jet  Results  with  HDL  Data 


Model  2  was  used  to  obtain  the  switch  times  of  a  laminar  attached 
jet.  The  supply  pressure  was  fixed  at  pj.  =  0.0333  kPa  so  that  the 
Reynolds  number  Re  =  1000  ,  based  on  nozzle  width  and  pj..  Tests  at 

1000  were  not  conducted  because  of  pressure  regulation  difficul- 
very  low  supply  pressures.  Based  on  the  work  of  Comparin  et 
al  and  the  quiet  hot-film  probe  output,  lamunar  flow  was  presumed. 

Low  level  control  inputs  were  also  difficult  to  regulate;  therefore, 
the  imposition  of  an  input  pressure  ramp  that  increased  to  a  final 
level  (as  used  in  the  turbulent  tests)  was  discarded.  Instead  the 
laminar  input  signal  was  adjusted  so  that  switching  occurred  on  a 
rising  ramp.  Because  the  laminar  control  signal  was  different  from  the 
turbulent  signal,  the  presentation  of  the  data  is  slightly  different. 

For  these  tests,  the  slope  of  the  ramp,  Ap/At,  is  plotted  against  sv/itch 
time.  Figure  27  shows  the  laminar  experimental  results  in  good  agree¬ 
ment  with  the  theory.  A  single  theoretical  curve  is  given  because  the 
difference  between  the  analytical  results  for  open  and  blocked  inactive 
controls  was  insignificant. 

For  the  turbulent  jet,  analytical  results  of  attachment  point  dis¬ 
tance,  control  flow,  and  inactive  output  flow  were  presented  as  a  func¬ 
tion  of  time  (fig.  26) .  Except  for  the  attachment  point  distance,  the 
sp  variables  are  shown  in  figure  28  for  Model  with  a  laminar  jet. 

S....oe  the  attachment  point  is  initially  past  the  vent,  the  distance  xj 
(the  distance  from  the  supply  exit  to  the  intersection  of  the  jet  cen¬ 
terline  and  the  geometric  center). ine)  is  shown  instead.  The  results 
shown  in  figure  28  are  for  a  ramp  input  APca/^t  =  7.4  x  10-“.  The  dis¬ 

tance  Xj  starts  from  a  negative  value  (jet  deflection  is  slightly  nega¬ 
tive)  and  increases  smoothly  until  it  approaches  the  splitter  (xj  =  10) . 
The  control  flow  Qca  increases  monotonically  after  an  initial  lag.  From 
a  practical  viewpoint,  the  output  flow.s  Qca  and  Qoi  are  of  interest. 

Note  that  there  is  always  some  flow  from  the  inactive  output  Qoi/  and 
unlike  the  turbulent  case,  there  is  not  as  large  a  change  in  either  out¬ 
put  flow  with  switching. 


4.5  Comparison  of  the  Switch  Time  of  a  Laminar  and  a  Turbulent  Jet 

Model  2  with  a  splitter  and  the  inactive  control  open  was  used  to 
obtain  laminar  and  turbulent  analytical  data.  Since  the  intermediate 
theory  has  been  shown  to  be  valid,  the  analytical  results  for  laminar 
and  turbulent  switch  times  can  be  compared.  Turbulent  values  are  taken 
from  the  data  used  in  figure  20,  but  are  replotted  as  Pcg  -  p^y  in  figure 
29.  The  theoretical  input-control  signals  are  also  shown.  For  this  com- 
the  values  of  laminar  data  of  (fig.  27)  are  used.  From  the  laminar  val¬ 
ues  of  the  slope,  Apramp/At,  corresponding  to  the  slope  of  the  turbulent 
control  ramp  inputs,  the  laminar  switching  times  are  determined.  Com¬ 
parison  of  the  laminar  and  turbulent  switching  curves  on  figure  29  shows 
that  the  normalized  laminar  sv.’itching  rate  is  faster  than  the  turbulent. 
Figure  29  also  shows  representative  real-time  values  for  switching. 
Laminar  real  times  range  from  15  to  22  ms,  whereas  the  turbulent  values 
range  from  3  to  18  ms. 

4 . 6  Response  of  a  Flueric  Wall-Attachment  Amplifier 
TMost  Complex  Modell 

By  allowj.ng  the  jet  to  reattach  to  the  opposite  wall  after  separat¬ 
ing  from  the  initial  wall  in  the  analytical  model,  the  complete  time 
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Figure  27,  Laminar  switching  times,  Model  2. 


Figure  28.  Analytical  laminar  output  variables  during  switching,  Model  2. 
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Figure  29,  Comparison  of  laminar  and  turbulent  switch  times. 


history  of  all  the  problem  variables  can  be  determined.  A  flueric  wall- 
attachment  amplifier  of  the  design  shown  in  figure  30  was  tested.  Im¬ 
portant  features  of  this  particular  amplifier  are:  (1)  the  control  re¬ 
sistances  do  not  depend  directly  on  the  spacing  between  the  jet  and  the 
wall  when  the  corner  is  cut  away  because  the  smallest  area  through  which 
the  control  flow  passes  into  the  bubble  is  the  control  nozzle;  (2)  the 
measured  output  pressure  is  the  dynamic  pressure  issuing  into  the  vented 
region;  (3)  the  offset  of  zero  and  the  wall  angle  is  0.21  rad;  and  (4) 
loading  the  outputs  has  no  effect  on  the  flow  regime  in  the  amplifier 
interaction  region  due  to  the  decoupling  action  of  the  large  vent  region. 

Figure  31  and  32  show  oscilloscope  traces  of  the  amplifier  response, 
and  figure  33  shows  the  comparison  between  the  predicted  and  the  meas¬ 
ured  response.  Note  that  the  predicted  response  does  not  ej:hibit  the 
overshoot  of  the  data.  This  overshoot  or  "ringing"  seems  to  be  attribu¬ 
table  to  the  Heliriholtz  response  of  the  output  passage  coupled  to  the 
transducer  cavity;  otherwise,  the  agreement  is  good.  It  is  important  to 
note  that  the  divergence  oj.  figure  19  does  not  occur  when  the  jet  attaches 
to  the  opposite  side.  This  shows  that  the  results  are  truly  indicative 
of  the  actual  mechanisms  occurring  in  an  amplifier. 
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Figure  30-  Fluerlc  bistable  wall  attachment  amplifier 


5 .  REMARKS  AND  CONCLUSIONS 

We  have  shown  that  the  two-dimensional  turbulent  and  laminar  jet 
analytical  models  agree  with  experimental  results  obtained  from  several 
different  flueric,  bistable  amplifier  geometries.  The  theory  allows  the 
calculation  of  the  steady-state  jet  attachment  point  location,  jet 
switching  times  with  and  without  a  splitter,  and  amplifier  output 
response. 

Application  of  the  present  theory  for  low  attachment-wall  offsets 
can  be  justified  because  of  the  agreement  with  Johnston's  data*"*  and 
the  results  presented  in  section  4.6,  Note  that  the  jet  centerline  in¬ 
tersected  the  opposite  wall  during  the  analytical  calculations  for  both 
cases.  As  noted  previously,  the  jet  touching  and  flowing  along  the 
opposite  wall  is  not  considered  as  a  separate  mode  of  switching.  This 
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Figure  33.  Comparison  of  theoretical  switching  response  with  experimentally  observed  data 


implies  that  attachment  to  the  opposite  wall  and  at  the  same  time  to 
the  original  attachment  wall  does  not  result  in  a  self-sustaining  pro¬ 
cess  ending  with  a  switch.  Muller's  data*®(fig.  10,  11,  and  14  of  ref¬ 
erence  15)  indicate  that  the  jet  will  return  to  its  original  position 
after  flowing  along  the  opposite  wall  unless  a  minimum  control  pulse 
duration  is  exceeded.  Keto‘‘  concluded  from  his  experimental  evidence 
that  the  attachment  bubble  must  be  broken,  or  vented  to  some  other  port 
before  a  switch  can  be  considered  self-sustaining.  Since  the  opposite 
wall  deforms  and  even  restrains  the  size  of  the  bubble  rather  than 
bursts  It,  the  present  model's  assumption,  where  the  bubble  shape  is 
altered  when  the  jet  flows  along  the  opposite  wall,  is  justified. 

The  agreement  obtained  for  the  laminar  jet  switching  is  particularly 
satisfying  as  the  present  trend  in  fluidics  is  for  low  power  consumption 
devices.  The  laminar  model  is  not  concerned  with  whether  or  not  the  jet 
is  actually  laminar,  but  rather  whether  or  not  the  analytical  laminar 
jet  velocity  profile  expression  is  suitable  for  use  in  predicting  the 
jet  response  at  the  lower  Reynolds  number.  The  laminar  numerical  re¬ 
sults  from  which  the  switching  characteristics  of  section  4.4  are  ob¬ 
tained  indicate  that  a  significant  portion  of  the  flow  is  directed  to¬ 
ward  the  opposite  output  line.  In  the  usual  bistable  device,  flow  out 
of  the  inactive  output  may  be  undesirable;  however,  if  the  amplifier  is 
operated  differentially,  it  may  be  of  little  consequence.  The  good 
laminar  results  indicate,  however,  that  laminar  devices  should  be  in¬ 
vestigated  and  that  the  present  analytical  model  can  be  used  to  describe 
them. 

‘‘Koto,  J.K.,  "Transient  Behavior  of  Bistable  Fluid  Elements,"  Proc.  2nd  Fluid 
Amplification  Symposium,  HDL,  Vol-  III,  May  1964. 


The  jet-spread  parameter  c  for  the  turbulent  model  is  the  only  pa¬ 
rameter  in  the  analysis  that  cannot  be  specified  arbitrarily.  The  value 
of  o  =  10  used  in  the  analysis  is  based  on  the  experimental  results  of 
others  and  is  as'.umed  constant  for  all  wall-attachment  amplifiers  (at 
least  for  the  range  of  aspect  ratios  of  one  to  3.3  of  the  models  for 
which  comparisons  have  been  made) .  The  general  agreement  of  the  theory 
indicates  that  the  use  of  one  constant  value  for  a  is  justified.  Chang¬ 
ing  0  as  much  as  20  percent  in  the  calculations  does  little  in  changing 
the  switching  time. 

From  the  designer's  viewpoint,  the  model  is  versatile  because  the 
switching  time  as  well  as  the  total  response  time  is  available  from  the 
computer  solution.  Analytical  results  for  all  other  variables  such  as 
flows  and  pressures  are  also  available.  From  computer  input  and  output 
data,  quantities  such  as  pressure  and  flow  recovery,  pressure  and  flow 
gain,  fan  out,  and  steady-state  characteristic  curves  are  obtainable  if 
desired,  and  the  present  results  show  that  this  is  feasible. 

Future  refinements  to  the  model  should  include  consideration  of  blunt 
and  cusped  splitter  fluid  amplifiers.  The  potential  of  the  present 
model  may  also  be  increased  by  including  asymmetrical  geometries,  curved 
attachment  walls  and  temperature  effects  through  known  relationships  be¬ 
tween  density,  viscosity  and  temperature. 


LIST  OF  SYMBOLS 


A  area  per  unit  depth 

AR  aspect  ratio  =  depth/nozzle 
width 

b  nozzle  width 

B  amplifier  half-width  at  the 
downstream  edge  of  the  control 

discharge  coefficient 
C  laminar  jet  coefficient 
d  distance  along  attachment  wall 
D  attachment  wall  offset  distance 
F  -force 
h  unit  depth 

k  jet  profile  coefficients 
K  resistive  parameter 
L  inductive  parameter 

m  vortex  velocity  function 
coefficient 

n  coordinate  perpendicular  to 
jet  centerline 
p  pressure 

Q  volumetric  flow  rate  per  unit 
depth 

r  radial  distance  associated  with 
bubble  vortex 

R  jet  radius  of  curvature 
Re  Reynolds  number  =  ~ 

s  coordinate  distance  measured 

along  jet  centerline  from  supply 
exit  plane 

t  time 

T  attachment  parameter 
u  velocity 

V  volume  per  unit  depth 
W  line  area  per  unit  depth 

X  coordinate  distance  along  axis 
of  unit 

y  coordinate  distance  normal  to  x 

Aye  displacement  of  jet  centerline  from 
geometric  axis  measured  at  down¬ 
stream  edge  of  control 

z  distance  along  a  streamtube 
a  attachment  wall  angle 
3  jet  deflection  angle 

angle  associated  with  vent 
restriction 

n  similarity  parameter-turbulent 
0  jet  attachment  angle 


\  subsidiary  attachment  variable 

V  kinematic  viscosity 

5  similarity  r^arameter-laminar 
p  density 

0  turbulent  jet  spread  parameter 
angular  coordinate 

(i)  angular  velocity  of  bubble 
vortex 

subscri pts 

a  active  side  (side  to  which  jet 
is  initially  attached) 

ap  attachment  point 
av  average  value 
b  attachment  bubble 
c  control 

ca  active  control  (control  on 
initial  attachment  side) 

ci  inactive  control  (control 
opposite  initial  attachment 
side) 

d  aynamic  (except  for  the  dis¬ 
charge  coefficient  C|j) 

e  entrainment 

em  associated  with  vortex  velocity 

i  inactive  side  (side  opposite 
initial  attachment  side) 

j  point  where  jet  centerline  in¬ 
tersects  unit  axis 

o  output 
oa  active  output 
oi  inactive  output 
q  actual  flow 
r  returned 
s  supply  ixit 
sp  leading  edge  of  splitter 
t  total  conditions 

V  vent 

va  active-side  vent 
vi  inactive-side  vent 
vo  virtual  origin  of  jet 

w  point  where  jet  centerline 
intersects  attaclimcnt  wall 

+  invariant  supply,  source,  or 
reference  condition 
superscript 

Prime  denotes  dimensional  form 
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APPFNOIX  A 


*  OSL/90  LISTING  FOR  MODEL  1,  MODEL  2  (LAMIN-'R  AND  TURBULENT), 

*  AND  MODEL  ). 

* 

SiJOS  I79302-T  ,  15, 50000, ORZEWl  RETURN  TO  BLDG  92 

4*  OSL/90  MOUNT  TAPE  S-72  (DSL90)  ON  "5,  SCRATCH  UN  A5 
VEXtCUYE  USER 

VDSL90 

tlEOIT  SYSCK2,SCHF4 

MBLDR  MAIN 

J.IEOIT 

* 

*  THIS  PROGRAM  SIMULATES  THE  EQUATIONS  GOVERNING  THE  TRANSIENT  RESPONSE 

♦  OF  A  BISTABLE  AMPLIFIER  WITH  OUTPUT  LINES 

« 

♦INITIAL  CONDITIONS  ARE  ARBITRARY-  FDR  THIS  PROGRAM  THETA0=0.6 


V0=6.0 

C)S&=1.0 

♦ 

QCl0=QC20=0VlO=0.O 

♦ 

QV2=. 2, 001=2.0 

« 

QO2=.5,PLEVEL0=.2 

♦THESE  CONDITIONS  WILL  GENERATE  THE  ACTUAL  INITIAL  CONDITIONS  FOR  THE 
♦PROBLEM  BY  T I HE= F I N T I M / 2. O 

♦  FINTIM  IS  TOTAL  RUN  TIME  IN  SECONDS.  FIRST  HALF  OF  FINTIM  IS  USED 

♦  TO  CALCULATE  THE  INITIAL  CONDITIONS  FOR  THE  PROBLEM  GIVEN,  THE 

♦  SECOND  HALF  IS  USED  TO  COMPUTE  THE  TRANSIENT  RESPONSE  TO  THE  GIVEN 

♦  INPUTS 

♦ 

♦  THERE  ARE  TEN  DIRECT  OUTPUTS 

♦  TIME! SECONDS) , DIMENSIONLESS  TIME,  REATTACHMENT  PO  I  NT , C ENTERL I NE  OF 

♦  JET  LOCATION, ACTIVE  CONTROL  FLOW , I NACT 1 VE  CONTROL  FLOW, ACTIVE  OUTPUT 

♦  LINE  FLOW, INACTIVE  OUTPUT  LINE  FLOW,ACTlVF.  OUTPUT  DYNAMIC  PRESSURE, 

♦  INACTIVE  OUTPUT  DYNAMIC  PRESSURE 

♦ 

♦BASIC  DIMENSIONAL  PARAMETERS 

♦  TIME’SECONDS 

♦  TR1SE=  RISE  TIME  OF  INPUT  PRESSURE  SIGNAL, SEC 

♦  TAUTR=THE  TRANSPORT  TIME  FOR  ONE  NOZZLE  WIDTH  BO  AT  PRESSURE  PO,  SEC 

♦  FLUID  CHARACTERISTICS 

♦  RHO=DENSITY,  K I LOGR AM/ CUB IC  METER 

♦  NU=KINEMATIC  VISCOSITY,  METER  SQRO/SEC 

♦  AMPLIFIER 

♦  SUPPLY 

♦  0O=SUPPLY  NOZZLE  WIDTH,  MEIERS 

♦  802=B0  SQRD,  MgTERS  SQRO 

♦  PO=SUPPLY  PRESSURE,  KILOPASCAL 

♦  UPLIIS  =  SUPPLY  REFERENCE  VELOCITY,  M/SEC 

♦  CCEF  44,72  =  C0NVeRS ION  FACTOR  UF  SORT ( 2 . ♦ 1 000 ) 

♦  CONV-REF  VOL  FLOW  PER  UNIT  DEPTH  FOR  0.  EXIT  PRES,  METER  SQRO/SEC 

♦  INPUTS 

♦  Tl,T2=TIME  CONTROL  INPUTS  ARE  APPLIED,  SECONDS 

♦  Z1,Z2®F0RCING  FUNCTION  UF  UNITY  AMPL  APPLIED  AT  TIME  Tl,f2 

♦ 

♦BASIC  N0N0IMEN910NAL  PARAMETERS 

♦  TIME 

♦  TMOND«TIHE  NONO I  MENS  I  ON AL 1 Z ED  WITH  TAUTR 
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*  JET  CHftRACTERISTICS 

*  RE  =  REYNOl.DS  NUMBER  BASED  ON  REFERENCE  VOLUME  FLOW 

*  REC1,REC2=  CONTROL  REYNOLDS  NUMBERS 

*  SIGMA^TURBULBNT  JET  SPREAD  PARAMETER 

*  C1,C2=LAMINAR  JET  COEFFICIENTS 

*  SO=UPSTREAM  DIST  TO  APPARENT  JET  SOURCE  LAMINAR  JET 

*  AMPLIFIER  GEOMETRY 

*  CLOSE! ,CLOSE2=0.  FDR  A  CONTROL  THAI  IS  OPEN  OR  HAS  A  SIGNAL 

*  CLOSE!, CLOSE2=l.  FOR  A  CONTROL  THAI  IS  CLOSED  AND  CAN  HAVE  NO  SIGNAL 

*  PRE!,.PRE2=0.  FOR  OPEN  CONTROLS  DURING  INCON  CALC 

*  PRE!,PRE2=!.  FOR  CLOSED  CONTROL  DURING  INCON  CALC 

*  POSTl ,POST2=0.  FOR  OPEN  CONTROL  DURING  TIME  RESPONSE  CALC 

*  POSTl,POST2=!.  FOR  CLOSED  CONTROL  DURING  TIME  RESPONSE  CALC 

*  ALPH=ATTACHM6NT  WALL  ANGLE  (RADIANS) 

*  D=SETaACK 

*  SC'CONTRDL  WIDTH 

*  BC2=CaNrROL  WIDTH  SQUARED 

*  XVI=DISTANCE  ALONG  ATTACHMENT  WALL  TO  VENT 

*  LGTHS=LENGTH  OF  SUPPLY  CHAMBER 

*  AREAS=iCRUSS  SECT.  AREA/UNIT  DEPTH  UF  SUPPLY  CHAMBER 

*  LGTHCsLENGTH  OF  CONTROL  CHANNEL 

*  AKEAC=CRaSS  SECT.  AREA/UNIT  DEPTH  OF  CONTROL  CHANNEL 

*  LGTHV=LENGTH  OF  VENT  CHANNEL 

*  AREAV=CRaSS  SECT.  AREA/UNIT  DEPTH  OF  VENT  CHANNEL 

*  LGTHR=  length  OF  OUTPUT  LINES 

*  AREAR=  CROSS-SFCTIONAL  AREA  OF  OUTPUT  LINE 

* 

OPROBLEH  VARIABLES 

*  VELOCITY 

*  UE=BU8BLE  VORTEX  DRIVING  VELOCITY 

*  UE2=S0UARE  OF  BUBBLE  VORTEX  DRIVING  VELOCITY 

*  UEH=iMAX  VORTBX  DRIVING  VEL  BASED  ON  ENTRAINMENT  STREAMLINE 

*  M!,M2=MATCHING  COIF  FOR  THE  VORTEX  VELOCITY  FUNCTION 

*  FLOW 

*  V=C0NTINUITY  EQ  INTEGRAL  FOR  THE  BUBBLE  VOLUME 

*  QS=SUPPLY  FLOW 

*  QS2=SUPPLV  FLOW  SQUARED 

*  QC!,QC2  =  CO''jrRL)L  FLOW 

*  OVl=VENT  FLOW 

*  «V2=  FLOW  THROUGH  THE  OPPOSITE  SIDE  VENT 
»  QRl=RETURNeD  FLOW 

*  QEIsENTRAINEO  FLOW 

*  U01Q02=  FLOWS  THROUGH  THE  OUTPUT  LINES 

*  OSO, VO, QC!0, QC20, 0V!0tUV20, 0010,0020  INIlIAL  CONDITIONS  FOR  INTEGRALS 

*  Q01AVL=  SUM  OF  ALL  FLOWS  ENTERING  AMPLIFIER  EXCEPT  BY  THE  ACTIVE 

*  OUTPUT 

*  PRESSURE 

*  PCI  ,PC2-C0NTR0L  PRESSUPxE  INPUT  FORCING  fUNCTIONS 

*  P1,P2=;PERCENT  OF  SUPPLY  PRESSURE  TU  CONTROL 

*  PIBI AS,P2BIAS=CCNTR0L  BIAS  LEVEL 

*  PB1=ATTACHED  SIDE  BUBBLE  PRESSURE 

*  PB2=UNATTACHeD  SIDE  PRESSURE 

*  DEI PB^PRESSURE  DIFFERENCE  ACROSS  JET=PB?-Pdl 

*  PAV=AVERAGS  PRESSURE  AT  EXIT  OF  SUPPLY  AND  CONTROLS 

*  FSY=  SUPPLY  NOZZLE  EXIT  PRESSURE 

*  P1B,P26=  CONTROL  EXIT  PRESSURES 
PV1  =  AMBIENT  PRESSURE  AT  VENT 

*  PV2=  PRESSURE  AT  OUTSIDE  OF  THE  OPPOSITE  SIDE  VENT 
Pai,P02=  EXIT  PRESSURES  OF  TKE  OUTPUT  LINES 

*  PR1SE=  TRANSIENT  VALUES  OF  INPUT  PRESSURE  DURING  THE  RISE  TIME 
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*  OPDT=  GRADIENT  OF  INPUT  PRESSURE  SIGNAL 

*  P0l,PD2=  DYNAMIC  PRESSURES  IMPINGING  ON  THE  OUTPUT  LINES 

*  POUri,POUT2=  TOTAL  PRESSURES  IMPINGING  3N  THE  OUTPUT  LINES 

*  POOUTl,PDOUT2=  DYNAMIC  PRESSURES  AT  THE  END  UF  THE  OUTPUT  LINES 

*  PLEV£L=  PRESSURE  LEVEL  IN  AMPLIFIER  MEASURED  AT  EDGE  OF  THE  BUBBLE 

*  VORTEX 

*  DP=  INCREMENTAL  PRESSURE  CHANGE  DUE  TO  “ISMATCH  BETWEEN  AVAILABLE 

*  FLCW  AND  demand  FLOW 

*  LINE  CHARACTERISTICS 

*  RCl,RC2=CONTROL  RESISTANCE 

*  RCMIN1,RCHIN2=  MINIMUM  CONTROL  RESISTANCES 

*  RCMAX=MAX  CONTROL  RESISTANCE  EASED  ON  ENTRAINED  WIDTH  OF  8  DEG  JET 

*  SPREAD  AT  CONTROL  EDGE 

*  RClOt RC2D=DIST  BETWEEN  JET  EDGE  AND  WAL.  I KC*DY: ) , ( KC-D YC )  AT  CONTROL 

*  01,02=  CONTROL  DISCHARGE  COEFFICIENTS 

*  012,022=  CONTROL  DISCHARGE  COEFFICIENTS  SQUARED 

*  CDMIN,CDMAX=  MINIMUM  AND  MAXIMUM  DISCHARGE  COEFFICIENTS  OVER  RANGE 

*  Or  OPERATION  OF  THE  CONTROLS 

*  CS=  SUPPLY  NOZZLE  DISCHARGE  COEFFICIENT 
«  RV1=VENT  RESISTANCE 

*  RV2=  RESISTANCE  OF  OPPOSITE  SIDE  VENT 

*  RQ1,RQ2=  RESISTANCES  OF  THE  OUTPUT  LINES 

*  LS=SUPPLY  inductance 

*  LC=CONTROL  INDUCTANCE 

*  LV=VENr  INDUCTANCE 

*  LV2=  inductance  OF  OPPOSITE  SIDE  VENT 

*  L0=  INDUCTANCE  OF  OUTPUT  LINE 

*  MOMENTUM 

*  FC1,FC2=  FORCES  ACTING  TO  DEFLECT  IHE  JET  AT  THE  CONTROLS 

*  FSY=  JET  MOMENTUM  FLUX  EQN  AT  THE  INTERACTION  REGION 

*  aETArOEFLECTION  ANGLE  DUE  TO  SUPPLY  AND  CONTROL  JET  INTERACTION 

*  ZeT=  ANGLE  WHICH  DETERMINES  POINT  AT  WHICH  MOMENTUM  IS  FIRST  PEELED 

*  OFF 

*  XSl=  ANGLE  DETERMINING  AMOUNT  OF  PEELED  OFF  MOMENTUM 

*  YS=  REMAINING  WIDTH  ON  OFF  SIDE  OF  JET  AFTER  PEELING  OFF 

*  TS=  MOMENTUM  REMAINING  UN  OFF  SIDE  OF  JET  TO  ARRIVE  AT  REATTACHMENT 

*  POINT 

T  =  NO  LOSS  MOMENTUM  EG  AT  ATTACHMENT  POINT 

*  B=  ARGUEMENT  OF  ROOT  OF  REATTACHMENT  POINT  EQN 

*  GEOMETRY 

*  R=RADIUS  OF  CURVATURE  OF  JET  CENTER  LINE  I  CL ) 

*  VCL=REATTACHMENT  bubble  VOL  ENCLOSED  BY  JET  CL 

*  XCL=DISTANCE  ALONG  WALL  TO  CENTER  LINE  ATTACHMENT  POINT 

*  XAP=DIST  ALONG  WALL  TO  RETURNED  FLOW  ATTACHMENT  POINT 

*  LSP=OIST  TO  INTERSECTING  OF  JET  CL  AND  AMPL  CL 

*  SPL=SPLITTER  DISTANCE 

*  COl==.DIFFEReNCE  BETWEEN  VENT  AND  ATTACHMENT  POINT  DISTANCE 

*  THETA=INTERSECTI0N  angle  OF  JET  CENTER  LINE  AT  ATTACHMENT 

*  Sl=ARC  LENGTH  OF  JET  FROM  EXIT  TO  ATTACHMENT 

*  SF=HALF  S1=1QC4T!NG  ARC  LENGTH  FOR  VORTEX  CENTER  LINE  IN  BUBBLE 

*  SEM=ARC  LENGTH  ASSOCIATED  WITH  THE  VELOCITY  UEM 

*  SV=ARC  LENGTH  OF  JET  CL  TO  RADIAL  LINE  THRU  VENT  EDGE 

*  ZETA=ANGLE  FOR  CALC  ARC  LENGTH  SV 

*  Y1=N0RMAL  OIST  FR  JET  CL  TO  RETURNED  FLOW  STREAMLINE  AT  ATTACHMENT 

*  YV=NORMAL  DIST  FR  JET  CL  TO  RETURNED  FLOW  STREAMLINE 

*  KC=DIST  FROM  EDGE  UF  CONTROL  TO  EDGE  OF  NO-SPREAD,  UNDEFLECTED  JET 

*  OYC=LATERAl  DISPLACEMENT  OF  JET  MEASURED  AT  CONTROL  EDGE 

*  DV=EFFECTIVE  RESTRICTION  FOR  VENT  FLOW 

*  X1,X2,=CL  COORD  WHERE  JET  RADIUS  INTERSECTS  OPPOSITE  WALL 

*  XIl=ANGLE  SUBTENDED  BY  SUPPLY  EXIT  POINT  AND  XI 
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*  XI2=ANGLE  SUeTENDED  BY  XI  AND  X2 
>*>  SWU  =  ARC  subtended  by  XI2 

*  LEW=OIST  ALONG  OPPOSITE  WALL  SUBTENDED  BY  XI2 

*  WL=RATIO  OF  LEW/SWO 

* 

*  COMPUTATIONAL  VARIABLES 

*  L=l.  LAMINAR  CASE 

*  L=0.  TURBULENT  CASE 

*  E=0.  IF  ATTACHMENT  POINT  HAS  NOT  REACHED  VENT  EDGE.  IF  IT  HAS,  E=l. 

*  RT=1.  TERMINATES  SIMULATION 

*  TTT  =  INITIALIiI.NG  CONSTANT 

*  TTTT=  INITIALIZING  FACTOR  FORINPUT  RAMP  Gr'aDIENT  CALCS 

*  IT=  ITERATIVE  COUNTER  IN  PRESSURE  LEVEL  CALCS 

*  tincon=time  the  incqns  are  calc 

*  SA,CA,CA2,TA»A1.A2,A3,AA,C3,C4,C7=CALCULATED  CONSTANTS 

*  GAM,SG,CG,ST,.CT,SB,CB,CB2=SU8SlOIARY  CALC.  OF  VARIABLE  TRIG  FUNCTIONS 

*  B2,BW,BW2,CW,.SUM1  ,CBR,C6L,SVI  =  SUBS1DIARY  ALGEBRAIC  AND  TRIG  CALC. 

*  SAD,GAD=ARGUMENT  AND  SUBSIDIARY  CALC  FOR  THETA  IMPLICIT  LOOP 

*  QVlI,ARGaVl,QVlDnT=  COMBINATION  FOR  CALCULATING  VENT  FLOW  OVl 

*  PLUS=  SIGN  OF  DIFFERENCE  BETWEEN  AVAILABLE  FLOW  AND  DEMAND  FLOW 

*  th£tao=initial  quess  for  theta  implicit  loop 

*  plevelq=  inital  guess  at  the  pressure  level  in  the  amplifier 

*  ARGN=SaRT  function  ARGUMENTS 

*  RQOTN=SORr  OF  ARGN 

* 

01401  FORMAT!/,  9H  INCON  V=,F8.4,7H  THE TA= , F8 . 4, 3H  R=,F8.4,5H  QC1=,F8.4, 
0  15H  0C2=,Fa.4,5H  QVl=,F8.4,4H  OS=,F8.4,/) 

NOSORT 

IFiriME.GT.O.)  GO  TO  10 
UPLUS  =  44. 72*SQRT<  PO/RHO) 

C0NV=60*UPLUS 

R£=C0NV/NU 

TAUTR=BO/UPLUS 

Pa2  =  l.-OSO>t>*2*(  I  .-1.  <'<2.*R0>  ) 

PAV«l.-QS0**2*( 1,+1./(2.*R0) ) 

PSY=,5*(PAV+PB2) 

802=60*80 
BC2  =  8C'^BC 
SA=SINI ALPH) 

CA=COS( ALPH) 

CA2=CA*CA 
TA=rAN( ALPH) 

Al=3. /SIGMA 
A2  =  (D«-.B)*CA 
A3  =  (D+.!j)*SA 
A4=i(0+.5)*TA 
C3=C1/C2 

C4=(C2/(2.*C1) )**3 
C7  =  .?)*C2/(2.*C  1)**2 
LC=2.*LGTHC/AREAC 
LS=2. OLGTHS/AREAS 
KC=D+BC*TA 

KCMAX  =  1,/S.141A*BC)*>!'2 

AREAR  =  A2<-SPL*S  A 

LO=2.*LGTHR/AREAR 

LV2=2.«LGTHV/AkEAV 

RDl=l./AREAR**2 

RD2=1 ./AREAR**2 

RV2=  1  ./ARf!AV**2 
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RVll=l./AftEAV**2 
IF(L.GT.O.)  GO  TO  9 
UEM-SQRT ( 1./3.  ) 

SEM=2./A1 
9  CONTINUE 

*■ 

^tPROGRAM  INITIAL  CONDITIONS 

TI-IO.OE5 

T2=;10.0E5 

TTT=0.0 

TTTT=0. 

TINCON=0.0 

CLQSEI=PREl 

CLOSE2=PR62 

WL*0. 

XI2=0. 

RT  =  0. 

PLEVEL=PL£VLO 

IT=iO.O 

m 

10  CONTINUE 

♦  START  CALC  OF  DISCHARGE  COEFF  AND  MIN  RESISTANCES 

* 

RECI=Aas (QCl+RE) 

REC2  =  ABS(GC2'*'RE) 

1F(REC1.LT.2500. )  GO  TO  I 

Dl=-(COMAX-CDMI.N)*R£Cl/UE44-l.25*COMAX-.25*CDMIN 
GO  TO  2 

1  IFIRECl.GT.lOOO. )  GO  TO  5 

5  Dl=CDHAX 

2  IF(ReC2.LT.2500.)  GO  TO  3 
02=-<C0HAX-CDMI,N)*KEC2/i.e4«-l.2  5*CDM«X-.25*C0MIN 
GO  TO  4 

3  IF(REC2.GT.1000. )  GO  TO  6 
GO  TO  4 

6  D2»C0MAX 

4  CONTINUE 
012--^01*0l 
022=02*02 

RCHINl=l ./(BC2*DIZ) 

RCMIN2=I./(BC2*022) 

♦ 

*  END  CALC  OF  DISCHARGE  COEFF  AND  MIN  RESISTANCES 

* 

*THE  FOLLOWING  GENERATES  THE  INITIAL  CONDITION  PRINTOUT  AND  STARTS 

♦  THE  TIRE  DEPENDENT  PROBLEM 

* 

IF  1TIME.lt. FINTIH/2. )  GO  TO  15 

IFITTT.GT.O. )  GO  TO  15 

TTT=l-0 

Tl*TIHE 

TZ^TIME 

TINCON=TIME 

CL0S£1»P09T1 

CL0S£2=P>.jlT2 

WRITE16,  I40U  V,THETA,R,QCl,0C2,QVl,0S 
15  CONTINUE 
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OEND  OF  TIME  LOOP 


*  START  CALC  OF  PRESSURE  LEVEL 

*■ 

IF.(TIME.L£.m  GO  TO  99 
00lAVL=QSFQCl+QC2  +  avi  +  «VlA4-0V2+Q02 
PLUS=SIGN(.l.  ,001AVL-q01 ) 

DP-Rni*p;.US*IQai-QOlAVL)**2 
PLEVEL=PLBVEL+OP 
IT*TIME 
99  CONTINUE 

* 

*  END  CALC  OF  PRESSURE  LEVEL 

* 

TNONO=( TIME-TINCQN)/TAUTR 
Z2=STEP( T2) 

m 

*  START  CALC  FOR  JET  DEFLECTION 

« 

QS2>QS«0S 

IF<QC2.GE.0.0)  FACTOR=l.O 

IFiQC2.LT. 0.0)  FACTQR^O-0 

FCl=( PAV+2.*(QC1/BC)**2) *BC 

FC2=( PB2+2  .*I OC2*FACTOR/BC)**2J*aC 

FSY=PSY-*-2.*QS2 

BETA=ATAN( {FC1-FC2)/FSY) 

* 

♦  END  CALC  for  jet  DEFLECTION 

4< 

GAM»( ALPHFBETA ) 

SG=»SIN(GAM3 

CG=COS(GAM) 

VCL»V^.5*R*R*( XI2-SIN(XI2) ) 

B2»A2*«2/(2.*VCL  +  A2>t-A4) 

* 

*  START  implicit  ROUTINE  FOR  ATTACHMENT  ANGLE  THETA 

* 

THETA»IMPL( THETAO, I.E-3.GAD) 

ST^^SINITHETA) 

CT=;COS<TH6TA) 

RaOTl  =  (BZ*(SG+.5*ST)-CG)**2-CG*CG+B2*(GAM  +  THETA>.5<‘SIN{2.*GAM5  ) 
803  SAD  =  ARCOS(-a2t'(SG  +  .5*ST)*CG-SQRT(  ABSiROOTl)  I  ) 

GAD=ABS( SAD) 

4i 

M  END  IMPLICIT  ROUTINE  FUR  ATTACHMENT  ANGLE  THETA 

« 

*  START  CALC  OF  GEOMETRIC  VARIABLES 

4l 

R=A2/ICG-CT) 

Sl=»iR*(GAM+;THETA) 

SB^RASINOETA) 

CB-R*COS (BETA! 

CB2=CB*CB 
IF(BETA)  29,29,20 
20  BW^CA2*(  A4«-TA*CB-SB) 
eH2=BW*aw 

CW=:A2*A2+2.A‘A2*CA*CB 

4c 

ASTATEMBNTS  2S-29  CALC  IF  JET  CL  INTERSECTS  OPPOSITE  WALL 


25  IF.(BETA.LT.ALPH)  GO  TO  29 
ARG5=BW2-CW 

IF(ARG5.LT.O.)  GO  TO  29 
«OOT5=SQRT( ARG5) 

X1=-BW-ROOT5 

X2=-BW*ROOT5 

IF(Xl.LT.O..OR.X2.LT.O.)  GO  TO  29 
XM=BETA--ARSIN(  {  SB-X1)/R) 

XI2=BETA-XI1+ARSIN( ( X2-SB1/R) 

SHO=R*XI 2 

LEW=2.*R*SINlXI2/2.) 

WL=L£W/SWO 

29  CONTINUE 

*  END  OPPOSITE  SIDE  INTERSECTION  EONS 

« 

*  start  calc  for  momentum  peeled  off  bv  splitter 

* 

ZET  =  ATAN(C8/(SPL“Sa)  I 
IF<ZET.GT,(PI/2.-TH£TA-ALPH))  go  to  31 
TS=l.O 
GO  TO  28 
31  CONTINUE 

SS=R*lBETA-ZET+;PI/2.  ) 

XSI  =  SQRT(C8->*2+.(  SPL-SB)**2) 

YS=?XSI-R 

* 

*  END  CALC  FOR  MOMENTUM  PEELED  OFF  BY  SPLITTER 

Ik 

*  START  CALC  OF  MOMENTUM  AT  ATTACHMENT  POINT 

* 

IFIL.GT.O. )  GO  TO  27 
S0=SIGMA/3, 

TS=sTANH(  (YS*SIGMA)/(SS+SOn 
GO  TO  28 

27  SO^C^fRE+QS 

TS*TANri(  YS*C2'»(gS*RE/(5SFSOn**{2./3.  ) ) 

28  CONTINUE 

B»l.5*(-TS<-(TS*ik3/3.  )  *2.  /  3.  ♦•(  2 .  /  3. +  75- (  TS**3/ 3.  )  )*COS(  THETA 
T»2.*C0Sl  (2.*PI-ARCOS(-BZ2.) )/3.  ) 

* 

•  END  CALC  OF  MOMENTUM  AT  ATTACHMENT  POINT 

« 

•  START  ENTRAINED  FLOW  AND  RETURNED  FLOW  CALCS 

• 

IF(L»  30f30.A0 

30  ARG6*AI*S1+1. 

911  R00T6«SURT ( ARG6) 

0Rl«,5fOS*R00T6*( l.-T) 

QE-l-^S^QS^IROOTB-U  ) 

GO  TO  49 
40  SUHI»S1>S0 
SUM1»S13-S0 

CttR«(QS2/REl**{ 1./3. ) 

0Rl»C3«iCBR*SUMl*k  (l./3.I*(i.-T) 

QEl”C3*C8R*SUMl*ik(  1./3.  )-.5*QS 
SEM“SO*( SORT! 27. )-l. ) 

UEM  =  Cl+  l  2.*C3-i./3.  I/SQRTO.IYQS 
49  CONTINUE 
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♦ 

* 

* 


ENO  ENTRAINED  FLDW  AND  RETURNED  FLOW  CALCS 


ARGa=CB2-8C2+2.*BC*Sfi 
917  ROOT8=SQRT( ARG8) 

DYC  =  -Ctt<  ROUTS 

4c 

♦BEGIN  CALC  FOR  CONTROL  RESISTANCE 

♦ 

RC10=1./ ( (KC+nYC)**2+l.E-5) 

RC2D=I./ < (KC-OYC) ♦♦2+l.E-S) 

rFIRClD.GE.RCMIN.AND.RClO.LE.RCMAX)  »C1=RC1D 
lE(RClD.GT.RCMAX)  RC1=RCMAX 
IFIRCIO.LT.KCMINI )  RCI=RCMIN1 

IF (RC2D.GE.RCM IN. AND.RC20.LE.RCMAX)  RC2=RC2D 
IF  (RCaU.GT.RCMAX)  RC2-=RCMAX 
IF(RC2D.LT.RCM!N2)  RC2=RCMIN2 
IF(RCI.EO.RCID)  Dl=a. 

IF(RC2.EU.RC2D)  D2=l. 

* 

♦END  CALC  FOR  CONTROL  RESISTANCE 

♦ 

♦BEGIN  CALC  FOR  VENT  RESISTANCE 

4c 

ZE.TA  =  ArAN(  (XVl*CA-S8)/(CS-D-.5-XWl  +  SA  )) 
SV=R*(8ETA+2ETA) 

1701  IFIL-GT.O.)  GO  TO  1704 

Y1  =  (A1^S14-1.  )/6.^ALQG<(l.+T)/(l.-TJ  ) 

YV=Y1*( Al^SV*l.)/(Al^Sl*l.) 

GO  TO  1708 

1704  Y1»C74( 1 .♦(Sl/SO+1. >)♦*( 2./3. I^ALOGI ( 1 . +T ) / ( I .-T ) ) 
YV  =  V1^(  <  SV  +  SO)  /  (  S1->S0  )  )^+{2./3.  ) 

1705  CONTINUE 
XCL=R+{SGfST)-A3 
XAP=XCL-Y1/SIN<THETAI 
C01=XAP-XV1 
E=CaMPAR(XAP,XVl) 

IFIE.GT.O.)  GO  TO  400 
DV-0. 

GO  TO  401 

400  DV=R-YV-( ( XV1^CA-SB)/SIN( ZETA) ) 

401  CONTINUE 
IFIE.GT.O.)  GO  TO  500 
RV1=0. 

GO  TO  502 

500  IFIC01.lt. AREAV)  GO  TO  501 
RVl=l./AREAV++2 

GO  TO  502 

501  RVl  -  1./ ( DV  +  +  2+1. E-5) 

502  CONTINUE 

♦ 

♦END  CALC  FOR  VENT  RESISTANCE 

«c 

♦BEGIN  CALC  FOR  BUBBLE  VORTEX  VELOCITY 

♦ 

SE-.5t‘Sl 

IFISE.LE.SEM)  GO  TO  602 
Ir-IL.GT.O.  )  GO  TO  700 

♦ 

♦VORTEX  DRIVING  VEL  BASED  UN  ENTRAINMENT  TURBULENT 
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*  * 


SEl  =  U/  t  Al*SE+l.  ) 

UE2  =  2.25*9ei  +  (l.-SEn*(l.-SEl  )*CJS2 
GO  TO  701 

m 

♦VORTEX  DRIVING  VEL  BASED  ON  CUBIC  EQ  LAMINAR  AND  TURBULENT 

♦ 

602  Ml=2.  +  (0S-UEMl/lSEH*t‘3) 

M2*-l.5*ML*SEM 

ue  =  Hl*SE**3  +  l'2*SE«*2+OS 

ue2=uE*ue 

GO  TO  701 

* 

♦VORTEX  DRIVING  VEL  BASED  ON  ENTRAINMENT  LAMINAR 

♦ 

700  CBL=(RE  +  Q97(SE+:S0)  (♦♦(I./3.) 

UE-Cl*OS+CaL^I I.-(C27{2.^C1 ) )*^2^tCBL*^2  )  ) 

Ue2=UE^U6 

701  CONTINUE 

♦ 

♦END  CALC  FOR  BUBBLE  VORTEX  VELOCITY 

♦ 

PB1=-.S^UE2+PLEVEL 
0ELPB=2.*0S2/R 
Pa2=PBl+DeLPB 
PAV  =  .5^(  PBUPB21 
PSY  =  .5*(  PAV<-PB21 

♦ 

♦  START  CALC  OF  INPUT  RAMP  GRADIENT 

♦ 

IF(TIME,LT.FINriM/2. )  GO  TO  16 
IF ( TTTT.GT.O. )  GO  TO  16 
DPDT=(P1-PAV) / (TR ISE/TAUTR) 

PB10=PAV 

TTTT=l. 

16  CONTINUE 

♦ 

♦  END  CALC  OF  INPUT  RAMP  GRADIENT 

♦ 

♦  control  RAMP 

« 

IF(TIME.GT.TRISE+FINTIM/2.)  GO  TO  7 
PRISE=(Pl-P01O)+( TIME-FINTIM/2. )/TRISE+PBlO 
GO  TO  B 
7  PRISE-PI 
8  CONTINUE 

IRtTIME.LG.FINTIM/2. )  PRISE=0.0 
PIB=PAV 
P28=P02 

START  CALC  OF  INPUT  AND  OUTPUT  PRESSURES 

PC1=( l.-CLOSEl )♦( PRISE+P18IAS)+CL0SE1*P1D 
PC2= ( 1.-CL0SE2 )♦( P2*Z2+P26I AS  INCLOSE 2^P2  0 
PDl=.75+0S2+<TStT-( T S+^3»T*^3 ) / 3. )/AREAK 
PD2=.  75X'gS2*(2,/3.-T5*TS  +  ^3/3.)/AREAR 
P0UT1=PDIFPB2 
P0UT2=P02FPB2 
PD0UT1==I  Q0I/AREAR)++2 
PD0UT2=I Qa2/AREAR)^^2 
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* 

*  END  CALC  OF  INPUT  AND  OUTPUT  PRESSURES 

* 

*  START  CALC  OF  INTEGRALS  FOR  VOLUME  AND  ALL  FLOWS 

« 

OS  =INTGRL(OSO  ,CQNV«( 1.0-PSY-OS*ABS10S)/CS*»21/(LS*B02)*CS) 
QCl-INTGRL(aC10«CONV*(PCl-PlB-RCl*t)C!*AtiStOCl  )  J  /  ( LC*0 02  )  *D1 ) 
QC2=INTGRL(aC20tC0NV*(PC2-P2B-RC2*QC2+AbS(0C2) ) / ( LC*8 0 2 ) *D2 ) 
LV=2.*lGTHV-*S0RT(  Kvn 
OVl  =  FCNSW(  E,0-  ,.0.  ,Qvin 

ARGevl  =  CONV+(  PVl-PBl-RVl^^QVl^ABSIOVl)  )/(LV*B02) 

QVIDOT=FCNSW(E,0. lO.tARGQVn 
OVLI=INTGRL(QV10,OV1DUT) 

V=INTGKL IVO.CONV* (QC1-0EL+QR1-KJV1)/B02) 

QV2^INrGRL{OV20,CONV*(PV2-PB2-RW2*OV2»ABSlQV2) )/tLV*B02) ) 
001=INTGRL(Q010,CONV*(PD1+Pfl2-P0l-R01*OU1*ABS(Q01) )/(L0*B02) ) 
QO2=lNTGRL<oa2O,CONV*(PO2-PO2-PB2-RO2*OO2*A0S( 002) ) /(L0*B02) ) 

♦ 

*  END  CALC  OF  INTEGRALS  FOR  VOLUME  AMD  ALL  FLOWS 

« 

*  SWITCH  CRITERION 

4) 

LSP=2-*SB 

IFIuSP.LT.SPL)  GO  TO  2000 
RT»l.O 

2000  CONTINUE 
SORT 

TITLE  RESPONSE  OF  A  TURBULENT  REATTACHED  JFT  IN  AN  AMPLIFIER  WITH  OUTPUT 

INTEiG  MILNE 

RELERR  Vn.OOl 

C0NTRLFINTIM=,04 

FINI.SH  ftT=l. 

CONST  PI.=  3.1416»  Cl*,4S43f  C2=.2752 

RANGE  CLOSEl,CLOSE2,0Cl.QC2,PBl,PB2,PAV,PCl,RE,TAUTR,AREAR,OS,OO2,WL,... 

PLEveL,XCL,THeTAtUPLUS,0VUtQVl,QV2,V,SIGMA,DPDT 
PRINT  l.E-4,TNaND.XAP,LSP,QCl,QC2*Q01,Q02,PDDUTl,PDOUT2 
INCON  THETAO=0.60000f  V0=6.0000,  QSO=l.OOOO.  OC10=O.OOOOf 

aC20--0,0000.  QV10  =  0. 0000, 0V20=0. 2000,  0010=2.0000, 

0020=0.5000,  PLEVELO=0.2000 

* 

*  MODEL  I 

* 

PARAM 


PARAM 
PARAM 
END 

« 

*  FOR  MODEL  1  INPUT  CONTROL  AMPLITUDES  VARIED  FROM  0.  TO  0.5  FOR  BOTH 

*  INACTIVE  CONTROL  OPEN  AND  BLOCKED 

« 

*  MODEL  2 


L=  0  .  , 

PRE1=1. 

0=0.905 

PO^IO.O 

P1=.00, 

P2>^.00, 

B0=i2.  iE 

Pi ai AS= 

LGTKC=1 

SPL=10. 

TRISE=1 

CDHAX=. 


S  I.GKA=10.0, 
PRE2=I.0, 


RH0=1.2059, 
POSTING. 0, 


NU=1.4864E-5, 

P0ST2=1.0, 


,CS=0.85  , 


BC=1.0,  ALPH=.2C944,  XV1=  8.640, 

0,0,P2B iAS  =  0-U,PV 1  =  0.0, PV2=0.0,PO 1=0.0,  P02  =  0.  0  , 
0.0,AREAC=2.0,LGTHS=15.0,AREAS=3.0,LGTHV=10., AREAV=1.95, 
0  ^LGTHR=27.75 
.56-1 

9000, CDMIN=. 6000 
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O'  »J1  ^ 


TITLE  RESPONSE  OF  A  LAMINAR  REATTACHED  JET  IN  AN  AMPLIFIER  WITH  OUTPUT 


PRINT  5.E-3,TN0N0,XAP,LSP,QC1  3C2 , 00 1 , 002 ♦ PDOUT I , P OOUT2 
CONTRLFI.NT  I  M=4.0E--1 

PARAM  L=l.,  RHQ=l.205'3,  VU=  1 . 48  6  4E- 5 ,  ...  1 

PREl=1.0,  PRE2  =  i.O,  POSIl^O-'',  PJST2  =  l.O,  ...  2 

D- O.Sf  ...3 

P0=i3. 33336-2. CS  =  . 8, 

BO=^2.0E-3^  BC=1.0,  ALPH=  .20944 ,  XV1=  9.000. 

LGTHC=IO-0, ARE AC= 2.0, LOT HS= 15.0, ARE AS=3.0,LGTHV= 10-0. AREA V= 2.0,...  9 

SPL=10.  10 

END 


* 


*  THE  LAMINAR  MODEL  WAS  RUN  FOR  BOTH  INACT  IVL  CONTROL  OPEN  AND  BLOCKED 

*  FOR  RAMP  PRESSURE  SIGNALS  TO  Pl=.5  WITH  RISE  TIMES  FROM  30MS  TO  400MS 

* 

TITLE  RESPONSE  OF  A  TURBULENT  REATTACHED  JET  [N  AN  AMPLIFIER  WITH  OUTPUT 
PRINT  l.E-4, TnGND.XAP, LSP.QCI , 002,001,002, PDUUT1,PDUUT2 
CONTRLFINTIM=4.0E-2 

PARAM  L-0.,  SIGMA=10.0,  RHO=1.2059,  NU= 1 . 4864E-5,  ...  1 

PO=7.50,CS=0.85 

END 

* 

*  THE  TURBULENT  MODEL  WAS  RUN  FOR  BOTH  OPEN  AND  BLOCKED  INACTIVE 

*  CONTROLS  FOR  Pl=0.0  TO  0.5 


*■  MODEL  3 

♦ 


PARAM  TRISE=l.E-3 
PARAM  LGTHC=2.00 
PARAM  Pl=.25 
END 

* 

*  MODEL  3  WAS  RUN  FOR  PI-.25 


STOP 


SREMOVE  SYSCKl 
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APPENDIX  B 


5  *  AUTOSKEM  IS  A  PROGRAM  TO  DRAW  ELECTRONIC  SCHEMATICS  HOWEVER  IT  IS  ALSO 

6  *  USEFUL  FOR  DRAWINGS 

7  * 

8  ♦  THE  FOLLOWING  LISTING  IS  A  MODIFICATION  OF  THE  ORIGINAL  IN  ORDER  TO 

9  *  SHOW  THE  ACTUAL  TEST  SUE  OF  THE  MODEL 
10  * 


11 

*  ATTACHED 

JET  TEST 

MODEL  2  REV 

D  JMG 

30  DEC  71 

12 

♦  Or.’ENSInNS  ARE  IN 

CENTIMETERS 

13 

♦  COMPUTER 

RECOGNIZES  NUMBERS  AS 

INCHES 

SO  A  SCALE  FACTOR  IS  NEEDED 

14 

♦  ADJUST 

DRAWING  SIZE  ANO  CONVERT  DIMENSIONS 

15 

♦  FINAL  DRAWING  IS  ACTUAL  SIZE. 

SO  SCALE 

FACTOR  IS  1/2. 

54 

16 

* 

17 

FACTOR 

0,3937 

18 

ORIGIN 

0. 

0. 

19 

LINE 

.0000 

.0000 

.0000 

.  1000 

20 

LINE 

.0000 

.0000 

.1000 

.  0000 

21 

ORIGIN 

13.5 

9.5 

22 

LINE 

-5.0000 

.0000 

-5.2500 

.  0000 

23 

LINE 

-5.0000 

-.2500 

-5.0000 

.2500 

24 

LINE 

-5.0000 

.0000 

7-7500 

.0000 

25 

LINE 

7.5000 

-.2500 

7.5000 

.2500 

26 

LINE 

.0000 

.0000 

.0000 

.0000 

27 

LINE 

.0000 

.1000 

-.1000 

.1000 

28 

ARC 

-.1000 

.3500 

.2500 

210.000 

270.000 

29 

ARC 

-.9660 

-.  1500 

.7500 

30.000 

90.000 

30 

LINE 

-.9660 

.6000 

-4,0000 

.6000 

31 

CIRCLE 

-4.0000 

,00000 

o6000 

32 

LINE 

-4.0000 

-.6000 

-.9660 

-.6000 

33 

ARC 

-.9660 

,1500 

.7500 

270.000 

330.000 

34 

ARC 

1000 

-. 3500 

.2500 

90.000 

150. OOC 

35 

LINE 

-.1000 

-.1000 

.0000 

-. 1000 

36 

LINE 

.0000 

-.1000 

.0000 

-. 3425 

37 

ARC 

-.1000 

-.3425 

.1000 

300.000 

360.000 

38 

ARC 

.0000 

-.5157 

.1000 

120.000 

180.000 

39 

LINE 

-.1000 

-.5157 

-.1000 

-2.5157 

40 

CIRCLE 

.  1000 

-2.5157 

.2000 

41 

L  INE 

.3000 

-2.5157 

.3000 

-.5157 

42 

ARC 

.2000 

-.5157 

.1000 

.  000 

60.000 

43 

ARC 

.3000 

-.3425 

.1000 

180.000 

240.000 

44 

LINE 

.2000 

-.3425 

-2000 

-. 2426 

45 

LINE 

.2000 

-.2425 

1-7607 

-.5742 

46 

LINE 

1 . 7607 

-.5742 

1.7607 

-2.6157 

47 

L  INE 

1.7607 

-2.6157 

-0000 

-7.6200 

48 

LINE 

.0000 

-7.6200 

4.0000 

-7.6200 

49 

L  INE 

4.0000 

-7.6200 

2.1607 

-2.6157 

50 

LINE 

2.1607 

-2.6157 

2.1607 

-.6593 

51 

LINE 

2. 1607 

-.6593 

7.6200 

-1. 8197 

52 

LINE 

7.6200 

1.6197 

2.1607 

.6593 

53 

L  INE 

2.1607 

.6593 

2.1607 

2.6157 

54 

LINE 

2.1607 

2.6157 

4.0000 

7.6200 

55 

L  INE 

4.0000 

7.6200 

.0000 

7.6200 

56 

LINE 

.0000 

7.6200 

1.7607 

2.6157 

57 

L  INE 

1.7607 

2.6157 

1,7607 

.5742 

58 

LINE 

1.7607 

.5742 

.2000 

.2425 

59 

LINE 

.2000 

.2425 

.2000 

.  3425 

60 

ARC 

.  3000 

.  3425 

.1000 

120.000 

180.000 

61 

ARC 

.2000 

.5157 

.1000 

300.000 

360.000 

62 

LINE 

.3000 

-5157 

.  3000 

2.5157 

63 

CIRCLE 

.  1000 

2.5157 

.2000 

63 


64  Lir^E 

-.1000 

2.5157 

-.1000  .5157 

65  ARC 

.0000 

.5157 

.1000  180.000  240.000 

66  ARC 

-.  1000 

.3425 

.1000  .000  60.000 

67  LINE 

.0000 

.3425 

.0000  .1000 

68  circle 

-5.0800 

-6.9850 

.2381 

69  CIRCLE 

6.3500 

-6.9050 

.2381 

70  CIRCLE 

6.3500 

6.9850 

.2381 

71  CIRCLE 

-5.0800 

6.9850 

.2381 

72  LINE 

-7.62 

.0 

-7.62  -7.62 

73  LINE 

-7.62 

-7.62 

7.62  -7.62 

74  LINE 

7.62 

-7.62 
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Figure  H-1 .  Computer  drawing  of  J!DL  tent  model. 


